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Zusammenfassung 



Wir erstellen ein Modell, urn Hochenergiereaktionen von Hadronen und Pho- 
tonen zu berechnen. Der Gluon-Austausch beschreibt die perturbative und das 
Modell des Stochastischen Vakuums die nichtperturbative Wechselwirkung. Let- 
zteres fiihrt zum Quark- Confinement in Dipolen via eines Strings von Farbfeldern. 
Wir erforschen die QCD-Struktur der Dipol-Dipol Wechselwirkung im Impulsraum, 
vor allem die Wechselwirkung zwischen den Strings. Wir stellen den String als 
eine Ansammlung von stringlosen Dipoles dar, zeigen Confinement-Effekte in der 
Hochenergiestreuung und berechnen die unintegrierte Gluonverteilung von Hadro- 
nen und Photonen. Im StoBparameterraum des Streuprozesses untersuchen wir das 
Unitaritatslimit der S-matrix. Wir berechnen die StoBparameter-Profile der Hadron- 
Hadron und Photon-Hadron Streuung, bestimmen die Energiewerte bei welchen die 
Profile das Black-Disc-Limit erreichen, schatzen die stoBparameterabhangige Gluon- 
verteilung des Protons ab und diskutieren die Gluonsattigung. Wir vergleichen die 
Resultate des Modells fiir pp, irp, Kp, 7*p und 77 Reaktionen mit experimentellen 
Daten und identifizieren Sattigungseffekte in experimentellen MessgroBen. 



Abstract 



We develop a model to compute high-energy reactions of hadrons and pho- 
tons. The perturbative interaction is described by gluon exchange and the non- 
perturbative interaction by the stochastic vacuum model which leads to quark- 
confinement in dipoles via a string of color fields. We study the QCD structure of 
the dipole-dipole scattering in momentum space focussing especially on interactions 
between strings. We represent the string as a collection of stringless dipoles, show 
confinement effects in high-energy scattering and calculate unintegrated gluon dis- 
tributions of hadrons and photons. In the impact parameter space of the scattering 
process we investigate the unitarity limit of the 5-matrix. We calculate the impact 
parameter profiles for proton-proton and photon-proton scattering, determine the 
energy values at which the profiles saturate at the black disc limit, estimate the 
impact parameter dependent gluon distribution of the proton and discuss gluon sat- 
uration. We compare the model results for pp, irp, Kp, 7*p and 77 reactions with 
experimental data and identify saturation effects in experimental observables. 
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Chapter 1 
Introduction 



On the basis of the numerous tests in experiments, one can confidently say that 
Quantum Chromodynamics (QCD) is the correct theory of strong interactions. The 
asymptotic freedom of QCD has allowed us to understand the interaction of particles 
at high momentum transfers. The non-Abelian nature of QCD leading to color 
confinement prevents first-principle calculations of hadronic interactions in the non- 
perturbative region. In fact, it is a key issue in high-energy physics to understand 
and describe the non-perturbative structure of hadronic scattering processes. The 
unravelling of confinement effects in such interactions would be especially important. 

So far Lattice Gauge Theory [1] constitutes the only access to non-perturbative 
QCD physics from first principles. Numerical simulations of QCD on Euclidean 
lattices [2] give strong evidence not only for color confinement but also for chiral 
symmetry breaking and dynamical mass generation from the QCD Lagrangian. Un- 
fortunately, lattice QCD cannot be applied in Minkowskian space-time to simulate 
high-energy reactions since it is limited to the Euclidean formulation of QCD. 

Experiments show a rise of high-energy total cross sections with increasing center- 
of-mass (cm.) energy. The rise is slow for large-size particles (protons, pions, 
kaons, or real photons [3]) and becomes steep for small-size particles (highly virtual 
photons [4,5] or charmonia [6]) involved in the scattering process. It is highly 
unsatisfactory that we do not have a genuine understanding of the growing cross 
sections on the basis of the QCD Lagrangian. The most rigorous result we have is 
the Froissart bound [7] stating that hadronic cross sections cannot grow faster than 
ln 2 (s/so) for asymptotic cm. energies y/s. It is derived on the basis of very general 
principles such as unitarity and analyticity of the scattering matrix. Its derivation, 
however, does not provide any physical mechanism realizing the squared logarithmic 
increase. At present no quantum field theoretical understanding of the rising hadron 
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extension at high energies is available. 

It is hopelessly difficult to solve any of the mentioned problems of high-energy 
scattering from first principles alone. Consequently, models that approximate quan- 
tum chromodynamics are required. To be convincing, however, the proposed models 
should include the basic features of QCD, reproduce many experimental high-energy 
data and agree with lattice QCD simulations. 

In this work we develop a model to describe the high-energy scattering of hadrons 
and photons in the eikonal approximation [8]. Its central element is the gauge- 
invariant light-like Wegner- Wilson loop W r [C] with the boundary C and the rep- 
resentation r of SU(N C ) [1,9]. The scattering amplitude factorizes into the vac- 
uum expectation value of two correlated Wegner- Wilson loops (W ri [Ci]W r2 [C2]) and 
reaction-specific wave functions [10-13]. The Wegner- Wilson loops describe the path 
of color-dipoles and the loop-loop correlation function (W ri [Ci}W r2 [C2\) the dipole- 
dipole scattering. The most interesting are Wegner- Wilson loops in the fundamental 
and adjoint representation of £77 (3) which represent color-singlet quark- ant iquark 
dipoles and glueballs (adjoint dipoles), respectively. In our framework color-dipoles 
are given by the quark and antiquark in mesons or photons and in a simplified 
picture by a quark and diquark in baryons. The size and orientation of the color- 
dipoles in the hadrons and photons are determined by appropriate light-cone wave 
functions. In this sense the loop-loop correlation function constitutes the basis for 
a unified description of hadrons and photons. 

We evaluate the loop- loop correlation funtion (W ri [Ci]W r2 [C 2 ]) in the approach 
of Berger and Nachtmann [14] which resumes gluonic exchanges in the scattering 
process using a matrix cumulant expansion and the Gaussian approximation of the 
functional integrals in the gluon field strengths. The resumation is crucial to guar- 
antee the £-matrix unitarity and to investigate saturation effects in high-energy 
reactions. In additon to the interaction of two dipoles in the fundamental repre- 
sentation of £77(3) considered by Berger and Nachtmann [14], we compute also the 
interaction of a fundamental with an adjoint dipole in the representation of SU(N C ). 

We express the loop-loop correlation function (W ri [Ci]W r2 [C 2 ]) in terms of the 
bilocal gluon field strength correlators integrated over two connected surfaces. The 
surfaces enter by using the non-Abelian Stokes' theorem [11,15] to transform the 
line integrals over the gluon potentials in the Wegner- Wilson loops of the func- 
tional integral approach [10-13] into surface integrals over gluon field strengths. We 
use for the first time explicitly minimal surfaces, i.e., planar surfaces bounded by 
the Wegner- Wilson loops. In Euclidean space-time, this surface choice is usually 
used to obtain Wilson's area law [16]. We have recently shown that minimal sur- 
faces are actually required to ensure the consistency of our results for (W r [C]) and 
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(W^ ri [Ci]W ra [Cy) by using low-energy theorems [17]. The simplicity of the minimal 
surfaces is appealing: It allows us to show for the first time the analytic structure 
of non-perturbative interactions in momentum space. The previous pyramid mantle 
choice for the surfaces [13, 14, 18-23] did not allow such an analysis. 

We use in the gluon field strength correlator perturbative and non-perturbative 
correlations as needed to describe interactions of small and large-size particles. The 
perturbative correlator results from QCD while the non-perturbative correlator is 
modelled by the stochastic vacuum model (SVM) [16]. This combination allows us to 
describe short and long distance correlations in agreement with lattice calculations 
of the gluon field strength correlator [24,25]. Moreover, this two component ansatz 
leads to the static quark- ant iquark potential with color-Coulomb behavior for small 
and confining linear rise for large quark-antiquark separations as expected [26]. 

We use in the non-perturbative correlator the exponential correlation function 
adviced by lattice QCD investigations of long-distance correlations [25]. This cor- 
relation function stays positive for all Euclidean distances and, thus, is compatible 
with a spectral representation of the correlation function [27]. This represents a 
conceptual improvement since the correlation function that has been used in earlier 
applications of the SVM becomes negative at large distances [13, 14, 18, 20-23]. 

In spite of the numerous improvements our approach is still incomplete as it 
leads to energy-independent cross sections in contradiction to the experimental ob- 
servation. This is because of the missed gluon radiation in the present model. In 
this work, however, the energy dependence is introduced phenomenologically. In 
line with the experimental observation mentioned above, we ascribe a strong and 
weak energy dependence, respectively, to the perturbative and non-perturbative 
component. Assuming the same mechanism for the energy dependence of hadron 
and photon interactions (motivated by the successful two-pomeron fit of Donnachie 
and Landshoff [28,29]) we implement a powerlike energy dependence in the loop- 
loop correlation function (W ri [Ci]W r2 [C 2 ])- In this way a unified description of the 
energy behaviour of hadron-hadron, photon-hadron, and photon-photon reactions 
is obtained. The powerlike ansatz in combination with the multiple gluonic inter- 
actions resulting from the resumation method of Berger and Nachtmann is crucial 
to guarantee the Froissart bound [7]. The phenomenological energy dependence, of 
course, can only be an intermediate step. For a fundamental understanding of the 
energy dependence of cross sections quantum evolution has to be incorporated in 
our model. 

We adjust the model parameters to reproduce a wealth of high-energy scattering 
data, i.e. total, differential and elastic cross sections, structure functions and slope 
parameters for many different reactions over a large range of cm. energies. In 
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this way we have confidence in our model predictions for future experiments (LHC, 
THERA) and for energies beyond the experimentally accessible range. 

The non-perturbative component of our model gives confinement [16] due to flux- 
tube formation of color-electric fields between the color-sources in the dipole [17, 30, 
31]. The thickness of flux-tubes or confining QCD strings saturates for large dipole- 
sizes at the value of about one Fermi [17]. The general representation r of SU(N C ) 
kept for all computations leads to an exact Casimir scaling for static color-dipole 
potentials and QCD strings [17] in agreement with lattice QCD simulations [32,33]. 
In this work we show the intrinsic composition of QCD strings, the structure of 
string interactions, and manifestations of strings in high-energy scattering. 

We show for the first time the QCD structure of the perturbative and non- 
perturbative dipole-dipole interaction in momentum space within our model. We 
reproduce the known results for perturbative interactions between the dipoles, the 
two-gluon exchange [34,35], and give new insights into the non-perturbative dipole- 
dipole scattering process. This comes out as a sum of two parts: The first part 
describes the non-perturbative interaction between the quarks and antiquarks of 
the two dipoles and exhibits the same structure as the perturbative contribution. 
The second part represents the interaction between the strings of the two dipoles. 
The latter shows a new structure different from the perturbative two-gluon exchange. 

We find an extremely nice feature for the QCD string that confines the quark 
and antiquark in the dipole: The QCD string of length \f D \ can be exactly repre- 
sented as an integral over stringless dipoles of sizes £,\r D \ with < £ < 1 and dipole 
number density n(£) = l/£ 2 . This outstanding result is very similar to the perturba- 
tively computed wave function of a qq onium state in the large- N c limit where the 
numerous emitted gluons inside the onium state are considered as dipoles [36,37]. 

The decomposition of the QCD string into stringless dipoles allows us for the 
first time to extract the microscopic structure of the unintegrated gluon distribution 
of hadrons and photons J-h(x,k^) from our dipole-hadron and dipole-photon cross 
section via |fcj_| -factorization. We compare the unintegrated gluon distribution of 
the proton with those obtained in other approaches. 

The unintegrated gluon distribution of hadrons and photons J-~h ( basic, 
universal quantity convenient for the computation of many scattering observables 
at small x. It is crucial to describe processes in which transverse momenta are 
explicitly exposed such as dijet [38] or vector meson [39] production at HERA. 
Its explicit | k± \ - dependence is particularly suited to study the interplay between 
soft and hard physics. Moreover, the unintegrated gluon distribution is the central 
object in the BFKL [40] and CCFM [41] evolution equations. Upon integration over 
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the transverse gluon momentum \k±\ it leads to the conventional gluon distribution 
xGh(x,Q 2 ) used in the DGLAP evolution equation [42]. 

We show the manifestations of the QCD string explicitly in dipole-hadron cross 
sections at large dipole sizes and in unintegrated gluon distributions at small trans- 
verse momenta. We find further that the |fcj_|-factorization which is known in pertur- 
bative physics can be extended also to non-perturbative dipole-hadron interactions 
within our model. 

To study saturation effects that manifest the S-matrix unitarity, we investigate 
the scattering amplitude in impact parameter space where the S-matrix unitarity 
imposes the rigid black disc limit on the height of impact parameter profiles. We 
show explicitly that our profiles respect the black disc limit. Furthermore, the width 
of the impact parameter profiles is shown to increase logarithmically at asymptotic 
energies as needed to guarantee the Froissart bound [7]. 

We compute the impact parameter profiles for hadron-hadron and longitudinal 
photon-proton scattering for different cm. energies. We determine the energy values 
at which the impact parameter profiles saturate at the black disc limit for small 
impact parameters. The impact parameter profiles provide an intuitive geometrical 
picture for the energy dependence of the scattering process as they illustrate the 
evolution of the size and opacity of the interacting particles with increasing energy. 

We estimate the impact parameter dependent gluon distribution of the proton 
xG(x,Q 2 , \b±\) using its relation to the impact parameter profile for longitudinal 
photon-proton scattering. We find a \ow-x saturation of xG(x,Q 2 , \b±\) as a mani- 
festation of the S'-matrix unitarity. Consequently, the x-dependence of the integrated 
gluon distribution xG(x, Q 2 ) slows down from a powerlike to a squared logarithmic 
rise in agreement with complementary investigations of gluon saturation at low x. 

We compare the model results with experimental data and provide predictions for 
future experiments. We compute total cross sections <r tot , the structure function of 
the proton F 2 , slope parameters B, differential elastic cross sections da el /dt, elastic 
cross sections a el , and the ratios a el /a tot and a tot /B for proton-proton, pion-proton, 
kaon-proton, photon-proton, and photon-photon reactions involving real and virtual 
photons as well. The successful unified description of all these reactions indicates 
indeed the assumed universal pomeron contribution to the above reactions. 

Making use of the explicit saturation of the impact parameter profiles at the 
black disc limit, the energy values at which it is reached, and the logarithmic rise of 
the black disc radius, we show explicitly manifestations of S'-matrix unitarity limits 
in experimental observables. A squared logarithmic behavior of total cross sections 
of hadrons and photons is found to set in for cm. energies y/s > 10 6 GeV. The 
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ratios a el /a tot and a tot / B become universal for such energies, i.e., independent of 
the hadron species considered. For asymptotic energies, also the total hadronic cross 
sections become universal and increase in agreement with the Froissart bound [7]. 

The outline of the work is as follows: In chapter 2 we present the model and give 
the model parameters. In transverse momentum space considered in chapter 3 we 
show the structure of QCD string interactions, decompose the string into stringless 
dipoles, elaborate the total dipole-hadron cross section, and extract the unintegrated 
gluon distribution. We compare the latter with ones obtained from other approaches 
and use it also to compute the integrated gluon distribution of the proton. In impact 
parameter space investigated in chapter 4 we study the S'-matrix unitarity limits, 
calculate the impact parameter profiles for proton-proton and photon-proton scat- 
tering, determine the energy values at which the profiles saturate at the black disc 
limit, estimate the impact parameter dependent gluon distribution of the proton, 
and discuss gluon saturation. Finally, in chapter 5 we compare the model results 
with experimental data and identify saturation effects in experimental observables. 
The appendices contain conventions, hadron and photon wave functions, the ana- 
lytic continuation of the non-perturbative correlation functions from Euclidean to 
Minkowski space-time, and the non-forward scattering amplitude. 



6 



Chapter 2 

The Loop-Loop Correlation Model 



In this chapter we present the loop-loop correlation model. We describe briefly the 
functional integral approach to high-energy scattering in the eikonal approximation 
and compute in detail the expectation value of one Wegner- Wilson loop and the cor- 
relation of two Wegner- Wilson loops within a Gaussian approximation in the gluon 
field strengths. Perturbative dipole-dipole interactions (long-distance correlations) 
are described by perturbative gluon exchange and non-perturbative dipole-dipole 
interactions (short-distance correlations) are modelled by the stochastic vacuum 
model (SVM) [16]. We use the minimal surfaces to calculate the S"-matrix element. 
Finally, we introduce the energy dependence in a phenomenological way and specify 
the parameter values. 



2 . 1 Functional Integral Approach to High-Energy Scattering 



The T-matrix is the central quantity in scattering processes. It enters every ob- 
servable we intend to look at and is obtained from the S'-matrix by subtracting the 
trivial case in which the final state equals the initial state, 



where Pi and Pf represent the sum of incoming and outgoing momenta, respectively. 
Based on the functional integral approach to high-energy loop-loop scattering in the 
eikonal approximation [10-13], the T-matrix element for the reaction ab — > cd at 
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transverse momentum transfer (t = 



—q±) and cm. energy squared s reads 



x ^* c {z u ri) ip* d {z 2l r 2 ) 1- S rir2 (s,b x ,z 1 ,r 1 ,z 2 ,r 2 ) ip a (z 1 ,r 1 )i(j b (z 2 ,r 2 ) , (2.1.2) 



with the S rir2 matrix element 



W n [Ci]W n [C2] 

S rir2 (s,b±,z 1 ,r 1 ,z 2 ,r 2 ) = r — -, ^— . (2.1.3) 

W ri [Ci}) (W r2 [C 2 



G \ I G 

The crucial quantity in (2.1.3) is the light-like QCD Wegner-Wilson loop [1,9] 



W r [C] = Tr r Pexp 



-ig (j) dz»Ql(z)t a r 



(2.1.4) 



where the subscript r indicates a representation of SU(N C ), N c is the number of 
colors, Tr r = Tr r (- • ■) /Trll r is the normalized trace in the corresponding color-space 
with unit element ll r , g is the strong coupling, and Q^z) = G^(z)t® represents the 
gluon field with the SU(N C ) group generators in the corresponding representation, 
t", that demand the path ordering indicated by V on the closed path C in space- 
time. A distinguishing theoretical feature of the Wegner-Wilson loop is its invariance 
under local gauge transformations in color-space. Therefore, it is the basic object in 
lattice gauge theories [1,2,9] and has been considered as the fundamental building 
block for a gauge theory in terms of gauge invariant variables [43]. Physically, 
the Wegner-Wilson loop represents the phase factor acquired by a color-charge in 
the SU(N C ) representation r along the light-like trajectory C in the background 
field. In other words, the Wegner-Wilson loop describes a color-singlet dipole in the 
representation r of SU(N C ). 

Wegner-Wilson loops in the fundamental representation of SU (N c = 3) are espe- 
cially important since they represent color-singlet quark-antiquark dipoles. In this 
representation, the S rir2 matrix element (2.1.3) describes the elastic scattering of 
two quark-antiquark dipoles with transverse size and orientation r*j. The longitu- 
dinal momentum fraction of color-dipole i carried by the quark is Zj. The impact 
parameter between the dipoles is [19] 

b± = f lq + (1 - zi)r x - r 2q - (1 - z 2 )f 2 = r lcm - r 2cm , (2.1.5) 

where f iq (fig) is the transverse position of the quark (antiquark), r*j = 

and fi cm = Zifi q + (1 — Zi)Tiq is the center of light-cone momenta. In the eikonal 
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approximation to high-energy scattering the q and q paths form straight light-like 
trajectories. Figure 2.1 illustrates the (a) space-time and (b) transverse arrangement 
of the loops. 

The fi and distribution of the color-dipoles is given by the wave functions 
ipa t b and ipc^ that characterize the interacting particles. In this framework, the 
color-dipoles are given by the quark and antiquark in the meson or photon and in 
a simplified picture by a quark and diquark in the baryon. Concentrating in this 
work on reactions with a = c and b = d, only squared wave functions \ipi(zi, n)| 2 := 
ip*(z 1 ,f l )ifj a (z ll f 1 ) and \ip 2 (z2, r 2 ) | 2 := ^d( z 2, ^2) "06^2, ^2) are needed. We use for 
hadrons the phenomenological Gaussian wave function [23, 44] and for photons the 
perturbatively derived wave functions with running quark masses rrif(Q 2 ) to account 
for non-perturbative effects at low photon virtuality Q 2 [20] as discussed explicitly 
in Appendix B. 

The QCD vacuum expectation value (. . .)g in the S rir2 matrix element (2.1.3) 
represents functional integrals [11] in which the functional integration over the 
fermion fields has already been carried out as indicated by the subscript G. The 
model we use for the QCD vacuum works in the quenched approximation that does 
not allow string breaking through dynamical quark- antiquark production in color- 
dipoles in the fundamental representation of SU(3). The linear rise of static color- 
dipole potentials [16, 17] and dipole-hadron cross-sections at large dipole sizes as 
well as the l/|fc±| behavior of unintegrated gluon distributions at small transverse 
momenta \k±\ display explicitly the quenched approximation as shown in this work. 

In this work we perform computations with Wegner- Wilson loops in a general 
representation r of SU(N C ). This allows several investigations: (a) The large- N c 
limit of our scattering amplitudes can be examined which may give additional in- 
sights into the scattering process, (b) The high-energy scattering of color-dipoles in 
different representations describing different physical objects can be studied, (c) In 
Euclidean space-time, different representations of the Wegner- Wilson loop are ap- 
propriate to study the Casimir scaling hypothesis of static potentials or chromo-field 
distributions of color-dipoles [17,45]. (d) Since our model works in the quenched 
approximation, the investigation of the string breaking behavior of the static po- 
tential for color-dipoles in the fundamental and adjoint (glueballs) representation 
of SU(N C ) is possible: string breaking cannot occur in fundamental dipoles as dy- 
namical quark-antiquark production is excluded but should be present for adjoint 
dipoles because of gluonic vacuum polarization [17,46]. 
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Figure 2.1: High-energy dipole-dipole scattering in the eikonal approximation represented 
by Wegner- Wilson loops in the fundamental representation of SU(3): (a) space-time and 
(b) transverse arrangement of the Wegner- Wilson loops. The shaded areas represent the 
strings extending from the quark to the antiquark path in each color dipole. The thin 
tube allows to compare the field strengths in surface S\ with the field strengths in surface 
S2. The impact parameter b± connects the centers of light-cone momenta of the dipoles. 
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2 . 2 Vacuum Expectation Value of one Wegner- Wilson Loop 



To compute the expectation value of a Wegner- Wilson loop (2.1.4) in the QCD 
vacuum 



W r [C\) = (TirVexp 



-ig & dz»g a (z)t a r 



c 



G 



(2.2.1) 



we transform the line integral over the loop C into an integral over the surface S 
with OS = C by applying the non-Abelian Stokes' theorem [11, 15] 



W r [C}) G = (Tr r P 5 exp 



-i 9 - I da^{z)g; v {o,z;C zo ) t a r 



(2.2.2) 



G 



where Vs indicates surface ordering and o is an arbitrary reference point on the 



surface S. In Eq. (2.2.2), the gluon field strength tensor, G liu (z) 
parallel transported to the reference point o along the path C zo 

g^(o, z; C zo ) = $(o, z; C zo )- l g^(z)§(o, z; C zo ) 

with the QCD Schwinger string 



g%(z)t\ is 



<$>(o,z;C zo ) = Pexp 



g I dz^(z)t a r 

J Czo 



(2.2.3) 



(2.2.4) 



Due to the linearity of the functional integral, (Tr r . . .) = Tr r (. . .), Eq. (2.2.2) 
reduces to 



W r [C}) G = Ti r (P s exp 



■ 9 



i»- I da^(z)g; u (o 1 z;C zo )t ( l 



G 



(2.2.5) 



For the evaluation of (2.2.5), a matrix cumulant expansion is used as explained in [11] 
(cf. also [47]) 



Vs exp 



% | / da(z)g(o, z; C zo ) 
J s 



exp 



T- 



n=l 



>S 

i 9 -r 

n\ ' 2' 



G 



da(x 1 ) ■ ■ -da(x n ) K n (xx, 



(2.2.6) 



where space-time indices are suppressed to lighten notation. The cumulants K n 
consist of expectation values of ordered products of the non-commuting matrices 
g(o, z; C zo ). The leading matrix cumulants are 



Kx(x) = (g(o,x;C x )) G , 
K 2 (x 1 ,x 2 ) = (Vs[Q(o,x 1 ;C Xl )g(o,x 2 ;C X2 )]) G 

~2 ((G(o,x 1 ;C Xl )) G (g(o,x 2 ;C X2 )) G + (1 



(2.2.7) 
(2.2.8) 
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Since the vacuum does not prefer a specific color direction, K\ vanishes and K 2 
becomes 

K 2 (x u x 2 ) = (V s [S{o, zi; C xl )Q{p, x 2 ; C X2 )]) G . (2.2.9) 

Now, we approximate the functional integral associated with the expectation values 
(. . .)g as a Gaussian integral in the gluon field strength. Consequently, the cumu- 
lants factorize into two-point field correlators such that all higher cumulants, K n 
with n > 2, vanish 1 and (W r [C])G can be expressed in terms of K 2 



W r [C]) G = r Ti r exp 
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da^{ Xl ) / da pa (x 2 ) 
s Js 



V s [GUo, a*; C X10 ) t a r Q b (o, x 2 - C X20 ) t b r 



G 



'2.2.10) 



Due to the color-neutrality of the vacuum, the gauge- invariant bilocal gluon field 
strength correlator contains a 5-function in color-space, 



/ 9 

\4vr2 [&^(°> Xl > C *i°)&h( > X2 > Cx *°)] ) G =: \ 5abF ^P^ x ^ x 2, o; C Xl0 , C X2Q ) 

(2.2.11) 

which makes the surface ordering Vs in (2.2.10) irrelevant. The quantity F^ uprj will 
be specified in Sec. 2.4. With (2.2.11) and the quadratic Casimir operator C 2 (r), 



t a r t a r =t 2 r = C 2 (r)l r , 



Eq. (2.2.10) reads 



W r [C] ) = Tr r exp [-C 2 (r) X ss A] = exp 

G 



C 2 (r) 



Xss 



(2.2.12) 



(2.2.13) 



where 



Xss ■-- 



-iHl I da^( Xl ) [ da pu {x 2 )F^ pa { Xl ,x 2 , o; C XlQ , C X20 ) . (2.2.14) 
4 Js Js 



Our ansatz for the tensor structure of F pupa in Minkowski space-time - see (2.4.1), 
(2.4.6), and (2.4.3) - leads to xss = for light-like loops, as explained in Sec. 2.5. 
Consequently, the vacuum expectation value of the Wegner- Wilson loop becomes 



W r [C] 



1 . 



G 



(2.2.15) 



The deviation of ( W r [C\) from one is just a measure of the strength of parton 

splitting processes, e.g., q —> q + G. Our result (2.2.15) is consistent with the 
quenched approximation and remains valid for arbitrary surface choices. 

^^We are going to use the cumulant expansion in the Gaussian approximation also for perturba- 
tive gluon exchange. Here certainly the higher cumulants are non-zero. 
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In Euclidean space-time one obtaines Xss ^ since the loops are not light-like 
and F^ v is different from the Minkowskian one. This allows one to compute the 
potential of static color-dipoles [16, 17] from the expectation value of one Wegner- 
Wilson loop in Euclidean space-time as discussed briefly at the end of this chapter. 



2.3 The Loop-Loop Correlation Function 



The computation of the loop-loop correlation function (W ri [Ci] W t . 2 [C 2 ])g starts also 
with the application of the the non-Abelian Stokes' theorem [11, 15] that is used to 
transform the line integrals over the loops Ci >2 into integrals over surfaces Si j2 with 
9 Si 2 = C\ 2 as in the previous section 



W ri [C x \W n [C 2 ] ) g ={ Tr ri V s exp 



-i 9 - da^{xi)Gl v {oi^C xlol )t a r , 



x Tr r2 Vs exp 



9 



da pa (x 2 )Q b (o 2 ,x 2 \ 

C X2 o 2 ) t r2 



S 2 



G 



(2-3.1] 



where o\ and 02 are the reference points on the surfaces Si and 5*2, respectively, that 
enter through the non-Abelian Stokes' theorem. In order to ensure gauge invariance 
in our model, the gluon field strengths associated with the loops must be compared 
at one reference point o. Therefore, we require the surfaces Si and Si to touch at a 
common reference point 01 = 02 = o. 

Following the Berger-Nachtmann approach [14], the product of the two traces, 
Tr n (- • •) Tr r2 (- • •), over SU(N C ) matrices in the r*i and r 2 representation, respec- 



tively, is interpreted as one trace Tr r 



Tr T . 



•)/Tr r . 1 ig )r2 (1L 



that 



acts in the tensor product space built from the T\ and r 2 representations 



W ri [Ci]W r2 [C 2 



G 



Tr 



®r2 



P s exp[-z| lda^(xi)g; u (o, x x] C Xl0 ) t a ri ] 

1 L^. 



Si 



ll ri ® P 5 exp[-z| da pa (x 2 )G b pa (o,x 2 ;C X20 )t 
J s 2 



G 



11,. 



(2.3.2) 



Using the identities 

exp(^) ® 11 r2 = exp(i° <g> 1L 2 ) (2.3.3) 
Tl ri ® exp (** ) = exp ( ll ri ® ^ ) (2.3.4) 

the tensor products can be shifted into the exponents. With matrix multiplication 
in the tensor product space 

(t* ® ii r2 )« ® x 2 ) = titl ® n r 



(*« ® n r2 )(n n ® *J 



t a rl ® 



(2.3.5) 
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and the vanishing commutator 



(2.3.6) 



the two exponentials in (2.3.2) commute and can be written as one exponential 



W[d] W[C 2 ] } g = ( Tr ri55r2 V s exp 



-i- da' lv {x)g iJiV {o, x] C x 



G 



(2.3.7) 



with the following gluon field strength tensor acting in the tensor product space 



G; v (o,x;C xo )(t« ri ® ll r ,) for x e s, 

1 ^(o,x;a o )(ll ri ® t») for x G S 2 



(2.3.8) 



In Eq. (2.3.7), the surface integrals over S\ and S 2 are written as one integral over the 
combined surface S = Si + S 2 so that the right-hand side (rhs) of (2.3.7) becomes 
very similar to the rhs of (2.2.2). This allows us to proceed analogously to the 
computation of (W r [C])G in the previous section: With the linearity of the functional 
integral, the matrix cumulant expansion, the color-neutrality of the vacuum, and 
the Gaussian approximation now in the color components of the gluon field strength 
tensor Q^u(o, x; C xo ), only the n = 2 term of the matrix cumulant expansion survives, 
which leads to 



W rx [C x ]W n [C 2 



Tr ri r2 exp 



G 



(2.3.9) 



da^( Xl ) / da pa (x 2 )(Vs[G^(o, Xl ; C xlo )G pa (o, x 2 ; 
s J s 



Note that the Gaussian approximation on the level of the color components of the 
gluon field strength tensor (component factorization) differs from the one on the level 
of the gluon field strength tensor (matrix factorization) used to compute (W r [C]) 
in the original version of the SVM [16]. Nevertheless, with the additional ordering 
rule [48] explained in detail in Sec. 2.4 of [49], a modified component factorization 
is obtained that leads to the same area law as the matrix factorization. 



Using definition (2.3.8) and relations (2.3.5), we now redivide the exponent 
in (2.3.10) into integrals of the ordinary parallel transported gluon field strengths 
over the separate surfaces S\ and S 2 
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W n [d] W T2 [C 2 ] ) g = Tr ri ^ r2 exp 



dar{ Xl )\ da pa {x 2 )v s (gzM,^c xlo )gl(o,x 2 -c X20 ))jt a n ®t b J 



'2.3.10) 



G 



s-> 



y da^( Xl ) da»°{x 2 )V s (g"(o,x l] C xlo )gl(o,x 2 -,C X20 )) (t a ri ® t b J 



G 



\l da^(x t )l da»°{x 2 )V s (g"(o,x l] C xlo )g b (o,x 2 ,C X20 )) (t a r / ri ® fl 



G 



•S'2 



do^(xi) da p °{x 2 )V s (g«(o,x i; C Xlo )gUo,x 2] C X20 j) (fl n ® C,t» ) 



Here the surface ordering "Ps is again irrelevant due to the color-neutrality of the 
vacuum (2.2.11), and (2.3.11) becomes 



W^W^]) =Tr ri0r2 exp 

G 



(2.3.11) 



with 



XSiSj 



i^- [ da^( Xl ) [ da fKJ {x 2 )F )xvpa {x l ,x 2) o ] C Xl0) C X20 ) . (2.3.12) 
4 is, is, 



The symmetries in the tensor structure of -F^ pcr in Minkowski space-time - see (2.4.1), 
(2.4.6), and (2.4.3) - lead to XSiSi = Xs 2 s 2 = f° r light-like loops as explained in 
Sec. 2.5, and also to XSiS 2 = Xs 2 Si ='■ X- 2 Our final Minkowskian result for general 
SU (N c ) representations T\ and r 2 becomes 



W ri [C l }W T2 [C 2 \) =Tr ri ^ r2 exp 

G 



ix(t a ri ® K 2 



'2.3.13) 



After specifying the representations r \ and r 2} the tensor product t ri( ^ r2 :— 

can be expressed as a sum of projection operators Pj with the property Pjt ri0r2 = 



Aj Pj where 



^ri(g)r 2 — Aj Pj with Aj 



Tr ri (g) r2 (Pj 
Tr ri( g) r2 (Pi / 



(2.3.14) 



2 In Euclidean space-time a different result is obtained since XSiSi 7^ as shown in [17]. 
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which corresponds to the decomposition of the tensor product space into irreducible 
representations. 

For two Wegner- Wilson- loops in the fundamental representation of SU(N C ), 7*1 = 
r 2 = N c , that describe the trajectories of two quark- ant iquark dipoles, the decom- 
position (2.3.14) becomes trivial 

f» ® =Hkzlp ^±lp (2.3.15) 
with the projection operators 

(Ps)(a 1 a 2 )(l3 1 f3 2 ) = ~X (<Wi ^a 2 /3 2 + <W 2 <Wl ) (2.3.16) 
(Pa)(aia 2 )(/3i/3 2 ) = 2 ^ a 2/3 2 ~~ <^ai/3 2 ^a 2 /3i ) (2.3.17) 

that decompose the direct product space of two fundamental SU (N c ) representations 
into the irreducible representations 



N c <g> N c = (N c + l)N c /2 © N C (N C - l)/2 . (2.3.18) 

Using Tr Nc(g)Nc Hn c ®n c = N% and the projector properties 

P* a = P s ,a, Tr Nc0Nc P s = (iV c + l)iV c /2, and Tr Nc ^ Nc P a = (N c - l)N c /2 , 

(2.3.19) 

we find a simple expression for the loop-loop correlation function with both loops 
in the fundamental SU(N C ) representation 



^N c [C 1 ]W Nc [C 2 ]) G = ^±^exp 



2N r X 



N c — 1 



iVc+l ' 

-x 



2N C 

(2.3.20) 

and recover, of course, for N c = 3 the result from [14]. 

For one Wegner- Wilson-loop in the fundamental and one in the adjoint repre- 
sentation of SU(N C ), 7*1 = N c and r 2 = Nj? — 1, that can be used to describe the 
scattering of a quark-antiquark dipole with a glueball, the decomposition (2.3.14) 
reads 

Nil 
*& e ® ^-i = -yPi+2 P 2--P 3 (2.3.21) 

with the projection operators 3 Pi, P 2 , and P3 that decompose the direct product 



3 Thc explicit form of the projection operators Pi, P2, and P3 can be found in [50] but note 
that we use the Gell-Mann (conventional) normalization of the gluons. The eigenvalues, A*, of 
the projection operators in (2.3.21) can be evaluated conveniently with the computer program 
"Colour" [51]. 
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space of one fundamental and one adjoint SU(N C ) representation into the the irre- 
ducible representations 

N C ©N 2 -1 = N c © In c (N c - l)(N c + 2) © In c (N c + 1)(N c -2) (2.3.22) 

which reduces for N c = 3 to the well-known SU (3) decomposition 

3©8 = 3©15©6. (2.3.23) 

With Tr Nc(g)N 2_ 1 H]sr c (g)N 2 -i = ^c(N^-l) and projector properties analogous to (2.3.19), 
we obtain the loop-loop correlation function for one loop in the fundamental and 
one in the adjoint representation of SU(N C ) 



iV c 2 -l 



exp 



G 

:N C ■ 

l Y 

. 2 



(2.3.24) 



N c + 2 



I — Y 

2 



N c -2 
2(N C -1] 



exp 



l — Y 

2 



Note that our results given in Eqs. (2.3.20) and (2.3.24) are quite general since 
they are obtained directly from the color-neutrality of the QCD vacuum and the 
Gaussian approximation of the functional integrals. For the explicit computation of 
X-function in Eqs. (2.3.20) and (2.3.24) one has to specify the gluon field strength 
correlator Fp Vpa and the surfaces Si^ with the restriction dSip = Ci,2 that appear 
in the Y-f unc tion (2.3.12). This we do in the following sections. Furthermore, we 
concentrate in the following on the dipole-dipole scattering (2.3.20) as the relevant 
case and postpone the dipole-glueball scattering (2.3.24) for future investigations. 



2.4 Perturbative and Non-Perturbative QCD Components 

We decompose the gauge-invariant bilocal gluon field strength correlator (2.2.11) 
into a perturbative (P) and non-perturbative (NP) component 

F = F p 4- F NP (2 4 1) 

where gives the perturbative physics (short-range correlations) described by 

perturbative gluon exchange and F^ the non-perturbative physics (long-range cor- 
relations) modelled by the stochastic vacuum model (SVM) [16]. This combination 
allows us to describe long and short distance correlations in agreement with lat- 
tice calculations of the Euclidean gluon field strength correlator [24,25]. Moreover, 
this two component ansatz leads to the static quark-antiquark potential with color- 
Coulomb behavior for small and confining linear rise for large source separations in 
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good agreement with lattice data as shown in our recent work [17]. Note that besides 
our two component ansatz an ongoing effort to reconcile the non-perturbative SVM 
with perturbative gluon exchange that has led to complementary methods [52-54]. 

We compute the perturbative gluon field strength correlator F^ from the gluon 
propagator in Feynman-'t Hooft gauge 



d 4 k 



(2tt) 4 k 2 - m 



(2.4.2) 



G 



where we introduce an effective gluon mass mo to limit the range of the perturbative 
interaction in the infrared (IR) region. In leading order in the strong coupling g, 
the perturbative correlator is gauge-invariant already without the parallel transport 
to a common reference point and depends only on the difference z := X\ — 22, 



[IVpCT 



• \ 9 2 1 
Z) = V 2 2 

= - g l 

7T 2 











dz v {Z ° 9w Zp9 ^ + d Zfl 



{zpgua - z a g vp ) D P (z 2 ) 



d A k 



kyk^g^p kykpg^u + k^kpg va k^k a g up D p {k ) 

(2.4.3) 



with the perturbative correlation function 



d 



D' P (k 2 ) :=— / d A zD P (z 2 )e 



.ikz 



k 2 — rriQ 



(2.4.4) 



We take into account radiative corrections in perturbative correlations by replac- 
ing the constant coupling g 2 with the running coupling 

487T 2 

g 2 (z ± ) = 4vra s (4) = jtz — r a#2Wa2 T ^ 2A ^ 

(33-2iV / )ln[(|^| 2 + M 2 )/A 2 QCD \ 

in the final step of the computation of the eikonal function \. Here the renormaliza- 
tion scale is provided by \z\\ that represents the spatial separation of the interacting 
dipoles in transverse space 4 . In (2.4.5) Nf denotes the number of dynamical quark 
flavors, which is set to Nf = in agreement with the quenched approximation, 
Aqqd = 0.25 GeV, and M 2 allows us to freeze g 2 for \z±\ — > 00. 

If the path connecting the points x\ and X2 is a straight line, the non-perturbative 
correlator F^f depends also only on the difference z := %\— x-i. Then, the most gen- 
eral form of the correlator in four-dimensional Minkowski space-time that respects 

4 Only transverse separations appear in the final expression of \ as explained in Sec. 2.5. 
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translational, Lorentz, and parity invariance reads [12, 13] 

■= (*) + fZ ( ? (z) (2.4.6) 

^W* W = g^ 2 - 1) \dz~ ( z °9w~ z p9^) + ~fcT \ z p9va-z a g vp ) JD^z) (2.4.7) 



6(7V ( 



^\ J ^27r) 4 6 (j^uko-g^p k u k p gp a -\- kpk p g VCJ kpkug^p^j D\{k 



Fp V pa (z) — _ ^ (9/ipg^cr QliaQvp) D(z ) 

= 3( ^ 2 ! 1} - / 70 ^(^) • ( 2 - 4 - 8 ) 

and was originally constructed in Euclidean space-time [16]. In all previous appli- 
cations of the SVM, this form depending only on x\ and x 2 has been used. New 
lattice results on the path dependence of the correlator [55] show a dominance of the 
shortest path. This result is effectively incorporated in the model since the straight 
paths dominate in the average over all paths. 

In (2.4.7) and (2.4.8), a is the correlation length, G 2 := (£zG%(0)G% u (0)) is 
the gluon condensate [56], k determines the relative weight of the two different 
tensor structures, D and D\ are the non-perturbative correlation functions in four 
dimensional Minkowski space-time, and 

D[(k 2 ) := f d^zD^/a 2 ) e lkz . (2.4.9) 

(IK J 

Euclidean correlation functions are accessible together with the Euclidean correla- 
tor in lattice QCD [24,25]. We adopt for our calculations the simple exponential 
correlation functions specified in four dimensional Euclidean space-time 

D E (Z 2 /a 2 ) = D?(Z 2 /a 2 ) = exp(-|Z|/a) , (2.4.10) 

that are motivated by lattice QCD measurements of the gluon field strength cor- 
relator Fp N 1 f pa (Z) [24,25]. This correlation function stays positive for all Euclidean 
distances Z and, thus, is compatible with a spectral representation of the corre- 
lation function [27]. This means a conceptual improvement as compared with the 
correlation function used in earlier applications of the SVMwhich becomes negative 
at large distances [13, 14, 18,20-23]. By analytic continuation of (2.4.10) we obtain 
the Minkowski correlation functions in (2.4.7) and (2.4.8) as shown in Appendix C. 

The perturbative and non-perturbative gluon field strength correlators can be 
analytically continued from Euclidean to Minkowski space-time by the substitution 
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8 up — > — g^ p and the analytic continuation of the Euclidean correlation functions to 
real time, D^(Z 2 ) — > D^ ! (z 2 ), [12, 13]. An alternative analytic continuation has been 
recently proposed by Meggiolaro [57]. In our resent work [17] we have generalized 
Meggiolaro's analytic continuation from parton-parton to gauge-invariant dipole- 
dipole scattering. In this approach one considers the correlation of two Wegner- 
Wilson loops tilted by an angle O with respect to each other in Euclidean space-time 
and obtaines - after Maggiolaro's prescription for the analytic continuation of the 
angle - the S'-matrix element for dipole-dipole scattering in Minkowski space- 
time [17]. Both analytic continuations give, of course, the same results and allow us 
to study the effect of confinement examined in Euclidean space-time on high-energy 
reactions computed in Minkowski space-time [8,17,58]. In chapter 3 we show the 
important features of the confining QCD string and its manifestation in high-energy 
scattering. 



The Static Color-Dipole Potential 

To illustrate the meaning of the different contributions F^J}a c) (z) and F^J}a(z) 
in the non-perturbative correlator (2.4.6), let us consider the potential of a static 
color-dipole which is obtained from the expectation value of one Wegner- Wilson 
loop 5 computed in Euclidean space-time (E) [1,59] 

V r (R) = -hm Ln(W r E [C}) . (2.4.11) 

1 — >oo 1 

The computation of (Wj^C]} follows the same line as the one for the Minkowskian 
case shown in Sec. 2.2: The only difference is the non- vanishing Euclidean Xss * n 
the final step of the computation which has to be computed now with the Euclidean 
correlator . For a loop in the fundamental representation, V r = Vn c describes 
the static quark- ant iquark potential and for a loop in the adjoint respresentation, 
V r = Vjv2_! gives the static potential of a gluino pair (adjoint dipole). 

Considered in Euclidean space-time, the perturbative correlator (2.4.3) leads to 
the non-confining color- Yukawa potential of a static color-dipole [17], 

V r P (R) = -^§^eM-rn 9 R], (2.4.12) 

where C^r) is the quadratic Casimir operator defined in (2.2.12) and r denotes 
the representation of SU(N C ). The perturbative contribution dominates the static 
potential at small dipole sizes R. In the limit m g — > 0, Eq. (2.4.12) reduces to the 
well-known color- Coulomb potential [60]. 

5 Here the subtraction of the self-energy of the color-sources is understood. 
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The tensor structure of F lhV p£ c ' (z) given in (2.4.7) is characteristic for Abelian 
gauge theories, coincides with the tensor structure of the perturbative correla- 
tor (2.4.3), and does not lead to confinement when considered in Euclidean space- 
time since it gives the static color-dipole potential [17], 

V r NP nc (R) = - C 2 (r) ^I'T R 2 Um , (2-4-13) 

3{N* - 1) 

which vanishes exponentially for large dipole sizes R. 

MP 

In contrast, the tensor structure of F^ v ^\z) given in (2.4.8) can only occur in 
non-Abelian gauge theories and Abelian gauge-theories with monopoles. Its Eu- 
clidean version leads to the following static color-dipole potential [17] 

V r NP C (R) = C 2 (r) R 2 J*dp (R - p)pKi[p/a] (2.4.14) 

when computed with the minimal surface, i.e., the planar surface bounded by the 
loop as indicated by the shaded area in Fig. 2.1. For large dipole sizes, R > 0.5 fm, 
Eq. (2.4.14) reduces to a linearly increasing static potential 

= a r R + const. , (2.4.15) 

_R>0.5fm 

which leads to confinement [16, 17]. Therefore, we call the tensor structure in (2.4.7) 
containing (1 — k) non-confining (nc) and the tensor structurein (2.4.8) containing 
k confining (c). 

Color-confinement [16] is realized in our model by the formation of a flux-tube 
of color-electric fields between the color-sources in the dipole [17,30,31]. We study 
the intrinsic structure of this flux-tube or confining QCD string and the interaction 
between QCD strings in chapter 3. 

The QCD string tension a r is given by the non-perturbative confining compo- 
nent (2.4.8) and reads for a color-dipole in the representation r of SU(N C ) [16, 17] 

3/^ poo 3 2 

o^GAr^-^fl d Z'D(Z\a>) = C 2 (r)*-^, (2.4.16) 

with the exponential correlation function (2.4.10) used in the final step. Since 
the string tension can be computed from first principles within lattice QCD [26], 
relation (2.4.16) puts an important constraint on the three fundamental param- 
eters of the non-perturbative QCD vacuum a, G 2 , and k. With the values for 
a, G 2 , and k given in Sec. 2.7, that are used throughout this work, one obtains 



V r NP C (R) 
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for the string tension of the SU(3) quark-antiquark potential (r = 3) a value of 
(T3 = 0.22 GeV 2 = 1.12GeV/fm which coincides with those obtained from hadron 
spectroscopy [61], Regge theory [62], and lattice QCD investigations [26]. 

The Casimir scaling of the static potential can be directly seen from (2.4.12)- 
(2.4.14). It emerges trivially in our approach as a consequence of the Gaussian 
approximation explained in Sec. 2.2. Since the Casimir scaling hypothesis of the 
static potential has been verified to high accuracy for SU(N C ) on the lattice [33], 
this result has been interpreted as a strong hint towards Gaussian dominance in the 
QCD vacuum and, thus, as evidence for a strong suppression of higher cumulant 
contributions [63,64]. In comparison to our model, the instanton model can neither 
describe Casimir scaling [64] nor the linear rise of the confining potential [65]. 

The static potential of an adjoint dipole differs from that of a fundamental dipole 
only in the eigenvalue of the corresponding quadratic Casimir operator: C 2 (r) = 
C 2 (Nl -1) = N C for adjoint and C 2 (r) = C 2 (N C ) = (iV c 2 - 1)/(2N C ) for fundamental 
dipoles. Our model working in the quenched approximation has a shortcoming at 
large dipole sizes: string breaking cannot occur in fundamental dipoles as dynamical 
quark-antiquark production is excluded but should be present for adjoint dipoles 
because of dynamical gluon production. From Eqs. (2.4.15) and (2.4.16) it is clear 
that string breaking is neither described for fundamental nor for adjoint dipoles in 
our model. Interestingly, even on the lattice there has been no striking evidence 
for adjoint quark screening in quenched QCD [46]. It is even conjectured that the 
Wegner- Wilson loop operator is not suited to studies of string breaking [66]. 



2.5 Evaluation of the %- Function with Minimal Surfaces 

For the computation of the x-function (2.2.14) 

X • X "T" Xnc ~~ Xc 

-i^-f da^(x 1 ) [ da^(x 2 )^ pa + FZ ( a c) +F^) , (2.5.1) 

J Si J S2 

one has to specify surfaces S\^ 2 with the restriction dS\^ 2 = C\^ 2 according to the non- 
Abelian Stokes' theorem. As illustrated in Fig. 2.1, we put the reference point o at 
the origin of the coordinate system and choose for S\^ 2 the minimal surfaces that are 
built from the areas spanned by the corresponding loops Ci t2 and the infinitesimally 
thin tube which connects the two surfaces Si and S 2 . The thin tube allows us to 
compare the field strengths in surface Si with the field strengths in surface S 2 . 
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Due to the truncation of the cumulant expansion or, equivalently, the Gaus- 
sian approximation, the non-perturbative confining component Xc depends on the 
specific surface choice. This is not the case for the perturbative and the non- 
perturbative non-confining component. A different and more complicated result 
for x^ P was obtained with the pyramid mantle choice for the surfaces in earlier 
applications of the SVM to high-energy scattering [13,14,18,20-23]. In this work 
we use the minimal surfaces because of the following reasons. Minimal surface are 
usually used to obtain Wilson's area law [16,17]. The minimal surfaces are also 
favored by other complementary approaches such as the strong coupling expansion 
in lattice QCD, where plaquettes cover the minimal surface, or large- iV c investi- 
gations, where the planar gluon diagrams dominate in the large- N c limit. Within 
bosonic string theory, our minimal surfaces represent the worldsheets of the rigid 
strings: Our model does not describe fluctuations or excitations of the string and 
thus cannot reproduce the Luscher term which has recently been confirmed with 
unprecedented precision by Luscher and Weisz [68]. 

Internal consistencies of our model favorize also the minimal surface: Since our 
results for the vacuum expectation value (VEV) of a rectangular Wegner- Wilson 
loop lead to a static quark- ant iquark potential that is in good agreement with lattice 
data [17], we are led to conclude that the choice of the minimal surface is required 
by the Gaussian approximation in the gluon field strengths. Furthermore, we have 
shown in our recent work [17] that the minimal surfaces are actually needed to 
ensure the consistency of our results for the VEV of one loop (W r [C]) and the loop- 
loop correlation function {W ri [Ci]W r2 [C 2 ]). Phenomenologically, in comparison with 
pyramid mantles, the description of the slope parameter B(s), the differential elastic 
cross section da el /dt(s,t), and the elastic cross section a el (s) can be improved with 
minimal surfaces as shown in chapter 5. 

The simplicity of the minimal surfaces is appealing: It allows us to show for the 
first time the structure of non-perturbative dipole-dipole interactions in momentum 
space, to represent the QCD string as a collection of stringless dipoles and to com- 
pute unintegrated gluon distributions of hadrons and photons. This is shown in 
chapter 3 and [58]. 

In applications of the model to high-energy scattering [8,58,67] the minimal 
surfaces are the most natural choice to examine the scattering of two rigid strings 
without any fluctuations or excitations. Our model does unfortunately not choose 
the surface dynamically and, thus, cannot describe string flips between two non- 
perturbative color- dipoles. This could generate also the energy dependence in non- 
perturbative interactions. Recently, new developments towards a dynamical surface 
choice and a theory for the dynamics of the confining strings have been reported [54]. 
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Parametrization of the Minimal Surfaces 



The minimal surfaces Si and S 2 shown as shaded areas in Fig. 2.1 can be parametrized 
with the upper (lower) subscripts and signs referring to S\ (S2) as follows 



Si = { (x?(u,v) 

(2) I V (2) 



(< +un» +vr?),ue [-T,T], ve [0,1] 1 , (2.5.2) 

V(2 9 ) (e) (2)/ J 
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(2.5.3) 
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The infinitesimally thin tube is neglected since it does not contribute to the x~ 
function. The computation of the ^-function requires only the parametrized parts 
of the minimal surfaces (2.5.2), the corresponding infinitesimal surface elements 



da^ 



dx p dx u <9x M dx L 



du dv dv du 



du dv 



11 v u v 

n% n - rf 

(e) (2) (2) (e) 



du dv 



(2.5.4) 



and the limit T — > 00 which is appropriate since the correlation length a is much 
smaller (see Sec. 2.7) than the longitudinal extension of the loops. 



X C omput at ion 



Starting with the confining component 



xT 



% ^ / do^ixi) [ da pa {x 2 ) F^\z = Xl - x 2 ) 
4 JSi Js 2 
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da^fa) / da f,a (x 2 )(g w g ua -g IMa g up )tD(z 2 /a 2 ) , (2.5.5) 
Si Js 2 



one exploits the anti-symmetry of the surface elements, da pu = —dcr Ufl , and applies 
the surface parametrization (2.5.2) with the corresponding surface elements (2.5.4) 
to obtain 



NP 



X 



where 



2 (r*i • f 2 ) / dv\ I dv 2 lim / du\ \ du 2 iD(z 2 /a 2 ) , (2.5.6) 
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(2.5.7) 
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-ikz 



and the identities n e ■ r 2 = r\ ■ n e = and n ffi ■ n e = 2, evident from (2.5.3), have 
been used. Next, one Fourier transforms the correlation function and performs the 
U\ and u 2 integrations in the limit T — > oo 

lim / diti / du 2 iD(z 2 /a 2 ) 
T^ooJ__ T J_ T 

= I w^ iD{k2) ^ Jj ui i T du2e ~ 

= *D{k 2 ) expi-ik^ - r% + v x r{ - v 2 r%)] 6(k° - k 3 ) 5(k° + k 3 ) 

= - iD {2) (fig + uift - r 2q - v 2 f 2 ) , (2.5.8) 

where iD^ is the confining correlation function in the two-dimensional transverse 
space (cf. Appendix C) 

^ (2)(4) = / ^^^^ ■ (2 ' 5 - 9) 

The contributions along the light-cone coordinates have been integrated out so that 
Xc F is completely determined by the transverse projection of the minimal surfaces. 
Inserting (2.5.8) into (2.5.6), one finally obtains 

x c P = KAT2 2 T\ ' ^ / dvi / d v 2iD {2) (nq + Vin-r^-v^) . (2.5.10) 

^U V c ~ Jo Jo 

With D^(k±) obtained from the exponential correlation function (2.4.10), cf. Ap- 
pendix C, we find 

iD (2) (z ± ) = 2ira 2 [1 + (\z±\/a)} exp(-|4|/a) (2.5.11) 
which is positive for all transverse distances. 
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Continuing with the computation of the non-confining component 
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we exploit again the anti-symmetry of both surface elements to obtain 
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2 (n • k) (r 2 ■ k) - (fi • f 2 ) (k u - r)(F + k 6 ) %D\{k l ) (2.5.13) 



with z as given in (2.5.7). Again the identities n e • r 2 = ri ■ n e = and n e • n e = 2 
have been used. Performing the u± and u 2 integrations in the limit T — > oo, one 
obtains — as in (2.5.8) — two (^-functions which allow us to carry out the integrations 
over k° and k 3 immediately. This leads to 



N p 7T 2 G 2 (1 - K 
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where iD^ is the non-confining correlation function in transverse space defined 
analogously to (2.5.9). The V\ and v 2 integrations are trivial and lead (cf. Fig. 2.1b) 
to 



NP 



7T 2 G 2 



3(iV c 2 



1 — k) 



iD'} 2) (fig - f 2q ) + iD[ (Z> (fig - f 2l 



/(2) 



ZD'/^ (fg - fg) - iD'^' (fig - f ,g ) . (2.5.15) 
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Using .D^ (/c^), derived from the exponential correlation function (2.4.10) in Ap- 
pendix C, we obtain 



iD[ [2 \z ± ) = 7i a A [3 + 3(|4|/a) + (|4|/a) 2 ] exp(-|4|/a) 



(2.5.16) 



The non-perturbative components, \c P an d Xnc > ^ eac ^ to co ^ or transparency 
for small dipoles, i.e. a dipole-dipole cross section with <7D£>(fi, r 2 ) oc |fi| 2 |f 2 | 2 for 
|fi,2| - > 0, as known for the perturbative case [69]. This will be shown in Sec. 3.2. 
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The perturbative component \ P is defined as 
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k u k a g a p k p k p g^ a + k^k p g ucr k a k a g pp iDp(k 



and shows a structure identical to the one of Xnc given in (2.5.13). Accounting 
for the different prefactors and the different correlation function, the result for 
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NP 



2.5.15) can be used to obtain 
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where the running coupling g 2 {z±) is understood as given in (2.4.5). With (2.4.4) 
one obtains the perturbative correlation function in transverse space 
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2^ 



K (m G \z±\) 



(2.5.19) 



where Kq denotes the th modified Bessel function (McDonald function) 



In contrast to the confining component Xc F ' i the non-confining components, Xnc 
and x P i depend only on the transverse position between the quark and antiquark 
of the two dipoles and are therefore independent of the surface choice. 

Finally, we explain that the vanishing of XSxSi an d Xs 2 s 2 anticipated in Sec. 2.3 
results from the light-like loops and the tensor structures in F nupa . Concentrat- 
ing — without loss of generality — on Xs 1 s l i the appropriate infinitesimal surface 
elements (2.5.4) and the F^ p(T -ansatz given in (2.4.1), (2.4.6), and (2.4.3) are in- 
serted into (2.2.14). Having simplified the resulting expression by exploiting the 
anti-symmetry of the surface elements, one finds only terms proportional to nL, 



ri, and n e • z with z^ 



.1, I 



.1, ■-) 



(ui — u 2 )n^ + {yi — v 2 )r^. Since = and 



ri = 0, which is evident from (2.5.3), all terms vanish and XsiSi = is derived. 

Note that x — X P + Xnc + Xc P * s a real-valued function. Since, in addition, the 
wave functions \ipi(zi, fi)\ 2 used in this work (cf. Appendix B) are invariant under 
the replacement (fj — > — r», Zi — > 1 —Zi), the T-matrix element (2.1.2) upon inserting 
the results (2.2.15) and (2.3.20) becomes purely imaginary and reads for N c = 3 

T(s,t) = 2is fd 2 b ± e^ fd Zl d 2 ri fdz 2 d 2 r 2 \tjj 1 (z 1 ,f 1 )\ 2 \^ 2 {z 2 ,r 2 )\ 2 



1 2 
- cos - 
3 V3 



(2.5.20) 
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The real part averages out in the integration over r*j and Z{ since the x-function 
changes sign 



which can be seen directly from (2. 5. 10), (2. 5. 15) and (2.5.18) as {f\ — > —T\,z\ — > 
1 — Z\) implies f\ q — > r 15 . In physical terms, (r*j — > — 2, — > 1 — Zi) corresponds to 
charge conjugation i.e. the replacement of each parton with its antiparton and the 
associated reversal of the loop direction. 

Consequently, the T-matrix (2.5.20) describes only charge conjugation C — +1 
exchange. Since in our quenched approximation purely gluonic interactions are 
modelled, (2.5.20) describes only pomeron 6 but not reggeon exchange. 



2.6 Energy Dependence 

Up to now the T-matrix element (2.5.20) leads to energy-independent total cross 
sections in contradiction to the experimental observation. This is disappointing from 
the phenomenological point of view but not surprising since our formalism does not 
describe gluon radiation which would generate the energy dependence. Nevertheless, 
we introduce the energy dependence in a phenomenological way inspired by other 
successful models so that a unified description of hadron-hadron, photon-proton, 
and photon-photon reactions is achieved. The powerlike ansatz used for the energy 
dependence is crucial to guarantee the Froissart bound [7] as shown in chapter 4. 

Most models for high-energy scattering are constructed to describe either hadron- 
hadron or photon-hadron reactions. For example, Kopeliovich et al. [70] as well as 
Berger and Nachtmann [14] focus on hadron-hadron scattering. In contrast, Golec- 
Biernat and Wiisthoff [71,72] and Forshaw, Kerley, and Shaw [73] concentrate on 
photon-proton reactions. A model that describes the energy dependence in both 
hadron-hadron and photon-hadron reactions up to large photon virtualities is the 
two-pomeron model of Donnachie and Landshoff [28] . Based on Regge theory, they 
find hard pomeron trajectory with intercept 1 + thard ~ 1.4 that governs the strong 
energy dependence of 7*p reactions with high Q 2 and a soft pomeron trajectory 
with intercept 1 + e so f t ~ 1.08 that governs the weak energy dependence of hadron- 
hadron or 7*p reactions with low Q 2 . Similarly, we aim at a simultaneous description 

6 Oddcron C = —1 exchange is excluded in our model. It would survive in the following cases: 

(a) Wave functions are used that are not invariant under the transformation (f\ — > — fj, Zi — > 1 — z 2 ). 

(b) The proton is described as a system of three quarks with finite separations modelled by three 
loops with one common light-like line, (c) The Gaussian approximation that enforces the truncation 
of the cumulant expansion is relaxed and additional higher cumulants are taken into account. 
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of hadron-hadron, photon-proton, and photon-photon reactions involving real and 
virtual photons as well. 

In line with other two-component (hard + soft) models [21, 22, 28, 73, 74] and the 
different hadronization mechanisms in hard and soft collisions, our physical ansatz 
demands that the perturbative and non-perturbative contributions do not interfere. 
Therefore, we modify the cosine-summation in (2.5.20) allowing only even numbers 
of soft and hard correlations, (x P ) 2 ™ {x NP ) 2m with n, m G iV. Interference terms 
with odd numbers of hard and soft correlations are subtracted by the replacement 



cos [cy] = cos [c (x P + X NP )] -> cos [cx P ] cos [cx NP ] 



(2.6.1) 



where c = 1/3 or 2/3. This prescription leads to the following factorization of soft 
and hard physics in the T-matrix element (2.5.20), 



T(s,t) = 2is Jd 2 b ± e^ Jdz l d 2 r l Jdz 2 d 2 r 2 |^i(^i, fi)\ 2 \M^,r 2 ) 



2 (I p \ (I NP 
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3 \3 A / V3 A 



/ 2 p\ (2 NP 
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(2.6.2) 



In the limit of small x-functions, \x F \ <C 1 and Ix^l ^ 1) one S e ^ s 



T{s,t) = 2is d 2 b x e^ b± d Zl d 2 n dz 2 d 2 r 2 |^i(^i,ri)| 2 \^ 2 {z 2 , r 2 )\ 
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(x p ) 2 +(x NP ) 



NP\ 2 



(2.6.3) 



In this limit, the T-matrix element evidently becomes a sum of a perturbative and 
a non-perturbative component. We show in the following chapter that the pertur- 
bative correlations, {x P ) 2 , describe the well-known two-gluon exchange [34,35] to 
dipole-dipole scattering and the non-perturbative correlations, (x NP ) 2 y the corre- 
sponding non-perturbative two-point interactions. 

As the two-component structure of (2.6.3) reminds of the two-pomeron model of 
Donnachie and Landshoff [28], we adopt the powerlike energy increase and ascribe 
a strong energy dependence to the perturbative component x P an d a weak one to 
the non-perturbative component x NP 



P\2 



U p ) 



s R 



(X NP (S)) 2 := ( X NP y 



2 ^2 



s rf r 2 

S Rq 



(2.6.4) 
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with the scaling factor SqRq. The power like energy dependence with the exponents 



o ■ 

« e < e p < 1 guarantees Regge type behavior at moderately high energies, 
where the small- x limit (2.6.3) is appropriate. In (2.6.4), the energy variable s is 
scaled by the factor f 2 rf that allows to rewrite the energy dependence in photon- 
hadron scattering in terms of the appropriate Bjorken scaling variable x 

sf i 2 «4 = -, (2-6.5) 
Q z x 

where \fi\ is the transverse extension of the qq dipole in the photon. A similar factor 
has been used before in the dipole model of Forshaw, Kerley, and Shaw [73] and also 
in the model of Donnachie and Dosch [74] in order to respect the scaling properties 
observed in the structure function of the proton. 7 In the dipole-proton cross section 
of Golec-Biernat and Wiisthoff [71, 72], Bjorken x is used directly as energy variable 
which is important for the success of the model. In fact, also in our model, the f 2 r 2 2 
factor improves the description of 7*p reactions at large Q 2 . 

The powerlike Regge type energy dependence (2.6.4) can be derived in more the- 
oretical frameworks: A powerlike energy dependence is found for hadronic reactions 
by Kopeliovich et al. [70] and for hard photon-proton reactions from the BFKL 
equation [75]. However, these approaches need unitarization since their powerlike 
energy dependence will ultimately violate S'-matrix unitarity at asymptotic energies. 
In our model, we use the following T-matrix element as the basis for chapter 4 and 5 

T{s,t) = 2is [d 2 b x e^- $J - [d Zl d 2 ri [dz 2 d 2 r 2 |^i(zi,fi)| 2 |^ 2 ,r 2 )| 2 



(2.6.6) 



2 cos(^x P (s)) cos(lx NP (s)] - ^cos(^x P (s) ) cosf-j^l 



3 \3" v 7 V3" v 7 3 V3 /v v 7 V3 



Here the cosine functions in combination with the powerlike energy ansatz ensure the 
S'-matrix unitarity in impact parameter space and the Froissart bound [7] as shown 
in chapter 4. Furthermore, the multiple gluonic interactions associated with the 
higher order terms in the expansion of the cosine functions provide the mechanism 
which leads to saturation effects in cross sections at ultra-high energies. 

Having ascribed the energy dependence to the x-function, the energy behavior 
of hadron-hadron, photon-hadron, and photon-photon scattering results exclusively 
from the universal loop-loop correlation function Sdd. In this way a unified de- 
scription of hadronic and photonic interactions is achieved. 

7 In the model of Donnachie and Dosch [74], s \ f\ \ is used as the energy variable if both 
dipoles are small, which is in accordance with the choice of the typical BFKL energy scale but 
leads to discontinuities in the dipole-dipole cross section. In order to avoid such discontinuities, 
we use the energy variable (2.6.4) also for the scattering of two small dipoles. 
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2.7 Model Parameters 



Lattice QCD simulations provide important information and constraints on the 
model parameters. The fine tuning of the parameters was, however, directly per- 
formed on the high-energy scattering data for hadron-hadron, photon-hadron, and 
photon-photon reactions where an error (x 2 ) minimization was not feasible because 
of the non-trivial multi-dimensional integrals in the T-matrix element (2.6.6). 

The parameters a, k, G 2 , m G , M 2 , s Rq, e NP and e p determine the dipole- 
dipole scattering and are universal for all reactions described. In addition, there 
are reaction-dependent parameters associated with the wave functions which are 
provided in Appendix B. 

The perturbative component involves the gluon mass mg as IR regulator (or 
inverse "perturbative correlation length") and the parameter M 2 that freezes the 
running coupling (2.4.5) for large distance scales. We adopt the parameters 

m G = m p = 0.77 GeV and M 2 = 1.04GeV 2 . (2.7.1) 

The value of the gluon mass is important for the interplay between the perturbative 
and non-perturbative contribution. Using our perturbative component, it gives a 
reasonable "perturbative glueball mass" (GB) of Mq B = IvrtQ = 1.54 GeV. The 
value of the parameter M 2 in (2.7.1) that fixes the running coupling at large distances 
at a s = 0.4 is taken from [76]. 

The non-perturbative component involves the correlation length a, the gluon 
condensate G 2 , and the parameter k indicating the non-Abelian character of the 
correlator. With the simple exponential correlation functions specified in Euclidean 
space-time (2.4.10), we obtain the following values for the parameters of the non- 
perturbative correlator (2.4.6) 

a = 0.302 fm, k = 0.74, G 2 = 0.074 GeV 4 , (2.7.2) 

and, correspondingly, the string tension for the fundamental representation of SU(3) 
becomes 

(73 = - K = 0.22 GeV 2 = 1.12 GeV/fm , (2.7.3) 
18 

which is consistent with hadron spectroscopy [61], Regge theory [62], and lattice 
QCD investigations [26]. Using the above value for the correlation length, the non- 
perturbative component generates a "non-perturbative glueball mass" of Mq B = 
2/ a — 1.31 GeV that is smaller than the "perturbative glueball mass" Mq B and 
thus governs the long-range interactions as expected. 
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Lattice QCD computations of the gluon field strength correlator down to dis- 
tances of 0.4 fm have obtained the following values with the exponential correlation 
function (2.4.10) [25]: a = 0.219 fm, « = 0.746, G 2 = 0.173 GeV 4 . This value for k 
is in agreement with the one in (2.7.2), while the fit to high-energy scattering data 
clearly requires a larger value for a and a smaller value for G*2- 

The energy dependence of the model is associated with the energy exponents 
e NP and e p , and the scaling parameter sqRq 

e NP = 0.125, e p = 0.73, and s i$ = ( 47 GeVfm 2 f . (2.7.4) 

In comparison with the energy exponents of Donnachie and Landshoff [28, 29], e so f t ~ 
0.08 and 6h ar d ~ 0-4, our exponents are larger. However, the cosine functions in our 
T-matrix element (2.6.6) reduce the large exponents so that the energy dependence 
of the cross sections agrees with the experimental data as illustrated in chapter 5. 

We have fixed the above parameters by trying to reproduce as many experimental 
observables as possible. The numerous observables can, of course, only narrow the 
range of the parameters. Additional contraints have been lattice QCD simulations 
and other models. Nonetheless, since an error (x 2 ) minimization has not been 
possible, the given parameter values cannot be mandatory. 
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Chapter 3 



Anatomy of QCD Strings and 
Unintegrated Gluon Distributions 



In this chapter we investigate the QCD structure of the perturbative and non- 
perturbative dipole-dipole interaction in momentum space within our loop-loop cor- 
relation model. We reproduce the known results for perturbative interactions be- 
tween the dipoles and give new insights into non-perturbative interactions between 
QCD strings which lead to quark-confinement in dipoles. Non-perturbative string- 
string interactions show a new structure different from perturbative dipole-dipole 
interactions. The confining QCD string exhibits a very nice feature: It can be 
exactly represented as a collection of stringless quark-antiquark dipoles. This out- 
standing result allows us to extract the microscopic structure of unintegrated gluon 
distributions of hadrons and photons. Our unintegrated gluon distribution of the 
proton is compared with those obtained in other approaches. Confining QCD string 
manifestations are shown explictly in dipole-hadron cross sections and unintegrated 
gluon distributions. The (-factorization, known in perturbative physics, is found 
to be valid also in non-perturbative interactions within our model. 

The structural aspects exposed in this chapter have been possible because of the 
simple minimal surfaces (see Fig. 2.1) used in our model. The previous complicated 
pyramid mantle choice of the surfaces [13, 14, 18-23] did not allow such studies. 

We consider throughout this chapter only the lowest order contribution to dipole- 
dipole interactions, i.e., the small-x limit of the T-matrix element (2.6.3). With 
a phenomenological energy dependence assigned to the x-f unc ti° ns as discussed 
in Sec. 2.6, this limit can be used to describe experimental observables for x = 
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Q 2 /s > 10~ 4 . For lower Bjorken-x (higher cm. energy) values, the exact T-matrix 
element (2.6.6) is needed to ensure S-matrix unitarity conditions and to observe 
saturation effects in high-energy reactions as discussed in chapter 4 and 5. 



3.1 QCD Structure of Dipole-Dipole Scattering 
in Momentum Space 

The total dipole-dipole cross section in the small- x limit is obtained from Eq. (2.6.3) 
via the optical theorem 



s 



ImT(s ,t = 0) 



2 d 2 b± Idzx^rx dz 2 d 2 r 2 \'ip Dl (zi,r 1 )\ 2 \'ip D2 (z2,r2)\' 



X 9L 



NP i „,iVP\ 2 ' 



(3.1.1) 



where ^/so ~ 20 GeV denotes the cm. energy at which our model (2.6.2) reproduces 
the experimental observables. The dipoles in Eq. (3.1.1) have fixed Zi and \fi\ values 
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27T \fr. 



H\ri\ ~ \r Dl \) 5(zi - z Di ) 



(3.1.2) 



but are averaged over all orientations. Being interested in the momentum-space 
structure of the scattering, we use (2.5.9) to writte the ^-functions computed with 
minimal surfaces in the previous chapter, see (2.5.10), (2.5.15), (2.5.18), in the form 
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(2tt) = 



e ik ± (r lq -r 2q ) , e ik±(r lq -r 2q ) 



_ e ik ± (fi q -r 2q ) _ e ik±(r lq -r 2q ) 



(3.1.4) 
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(3.1.5) 
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with the Minkowskian correlation functions in transverse space 



ii)f\kl) = ^-Ur, (3.1.6) 

J_ G 

iD?\kl) = ^-t, (3.1.7) 

iD®(kl) = 12712 5 . (3.1.8) 

a (k 2 + a -2 ) 2 

The perturbative correlation function, Dp 2 \k±), results from Eq. (2.4.4) for k° = 
k 3 = 0, i.e., k 2 = —k 2 , which is enforced in the computation of x - functions. 
The derivation of the non-perturbative correlation functions, D'i (k±) and D^ 2 \kf), 
from Eq. (2.4.10) is shown in Appendix C. 

The component x P describes the perturbative interaction of the quark and anti- 
quark of one dipole with the quark and antiquark of the other dipole as evident from 
the f ig and dependence of (3.1.3) and Fig. 2.1b. The component Xnc nas the 
same structure as x P an d gives the non-perturbative interaction between the quarks 
and antiquarks of the dipoles. With the term quark used generically for quarks and 
antiquarks in the following, we refer to x P an d Xnc as quark- quark interactions. 

The component x^ P given in Eq. (3.1.5) shows a different structure. Here the in- 
tegrations over v\ and v 2 sum non-perturbative interactions between the gluon field 
strengths connecting the quark and antiquark in each of the two dipoles. These con- 
nections are manifestations of the strings that confine the corresponding quark and 
antiquark in the dipole and are visualized in Fig. 2.1b. Indeed, the non-perturbative 
component Xc P shows a flux tube (string) between a static quark- antiquark pair 
in Euclidean space-time [17,30,31] which confines the quark and antiquark in the 
dipole as explained in chapter 2. Therefore, we understand the confining component 
Xc P as a string-string interaction. 

The mixed contribution Xnc Xc P ^ na ^ occurs in the dipole-dipole cross sec- 
tion (3.1.1) gives the non-perturbative interaction of the quark and antiquark of 
one dipole with the string of the other dipole, i.e., it represents the quark-string 
interaction. 

To be able to examine the \k±\ - structure of dipole-dipole interactions, we carry 
out all integrations in Eq. (3.1.1) except, of course, the one over transverse mo- 
mentum that enters through the x-functions (3.1.3)-(3.1.5). The resulting per- 
turbative (P) and non-perturbative (NP) integrand of the total dipole-dipole cross 
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section 



d\k ± \ I y (s ,\k ± \) + I NP (s ,\k ± \) 



(3.1.9) 



d\k ± \ I F (soAk±\) + C(s ,\k±\) + lZ F (soA^\) ,(3.1.10) 



NPi 



show the following momentum-space structure: 

i p (s ,\k\) 

x 2 (^|1 - e^|^) 2 (^ 2 |1 - e«|^ fl2 ) 



^i-|^|(4vra s (^)) 2 [^(^^ 



(3.1.11) 
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(3.1.12) 



2 <^|1 - e^|^) 2 (^ 2 |1 - e^|V^) 
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24 



ifc? 



ib^\k\ 



(3.1.13) 



2 (VbJ tan 2 0!(l - e lk ^)\^ Dl ) 2 (^ 2 | tan 2 2 (l - e^ 2 )|^ 2 ) 
where the brackets (■ip Di \...\ip Di ) denote the averages 

<^|A# Di > = IdzdW^oX^n^Ai (3.1.14) 



and the dipole orientation fa is defined as the angle between transverse momentum 
k± and dipole vector f{. With Eq. (3.1.2) the integration over the dipole orienta- 
tions fa leads, respectively, to the Bessel function Jo(|*i| I^dJ) an d the generalized 
hypergeometric function 1 iF 2 (— 1/2; 1/2, 1; — fc 2 r 2 ./4) 



<^|l-e^|^> 



2 k 



-J q dfail-^ 

i-M\k\K\) , 



(3.1.15) 



(^|tan 2 ^(l-e^)l^) 



2vr J 



2tt 



dfatan 2 fa(l - e ik±r " D ') 



1 1 -kfr, 



2^2 



'1'Di 



(3.1.16) 



1 A review of generalized hypergeometric functions can be found in [77]. In the computer pro- 

'l 7 



gram Mathematica [78] "HypergeometricPFQ[{ — 1/2}, {1/2, 1}, — denotes this function 
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The important implications will be discussed in the following section. 

The integrand I p given in Eq. (3.1.11) describes the known perturbative two- 
gluon 2 exchange [34, 35] between the quarks and antiquarks of the two dipoles. The 
ingredients of I p are visualized for one combination of gluon exchanges in Fig. 3.1a: 
the paired horizontal lines represent the dipole factors (1 — e lkl - Tl ) that describe the 
phase difference between the quark and antiquark at separate transverse positions, 
the surrounding brackets indicate the average over the dipole orientations (3.1.14), 
the two curly lines illustrate the gluon propagator squared [iDp 2 ^] 2 , and the four 
vertices (dots) correspond to the strong coupling to the fourth power g A = (4ira s ) 2 . 

The integrand I NP = I^ p + I^ p given in Eqs. (3.1.12) and (3.1.13) describes 
the non-perturbative interactions: the quark-quark, string-string, and quark-string 
interactions identified by the appropriate correlation functions [if)'^} 2 , [iD^] 2 , 
and [iD[^ iD^]. These interactions are illustrated in Figs. 3.1b, 3.1c and 3. Id, 
respectively. Analogous to the perturbative interaction in Fig. 3.1a, the dashed and 
solid zig-zag lines represent, respectively, the non-confining (1 — re) ib± and the 
confining (n/k 2 ) iD^ 2 ' non-perturbative correlations, the shaded areas symbolize the 
strings, and the four vertices (squares) in each figure indicate the "non-perturbative 
coupling" to the fourth power g% P := (7r 2 G2/24) . 

The integrand I^ p describes the non-perturbative interactions between the quarks 

and antiquarks of the two dipoles and exhibits the same dipole factors (1 — e lkj - n ) 
that appear in the perturbative integrand (3.1.11). I^ p contains three components: 
the non-confining component, the confining component, and their interference term. 
While the non-confining component visualized in Fig. 3.1b has the same structure 
as the perturbative contribution, see also Eqs. (3.1.3) and (3.1.4), the confining 
component shown in its more general form in Fig. 3.1c comes from the interaction 
between the quarks and antiquarks at the endpoints of the strings, which will be 
further discussed below Eq. (3.1.17). The interference term describes the quark- 
string interaction as illustrated in Fig. 3. Id. Note that it reduces entirely to an 
interaction between the quarks and antiquarks of the dipoles with the additional 
denominator \jk\ generated by the integrations over the variables V\ and v-i in the 
confining component \c P given in Eq. (3.1.5). 

The integrand I^ p describes the non-perturbative string-string interaction shown 
in Fig. 3.1c. The new angular dependencies in the string-string interaction, the mod- 
ified dipole factors tan 2 ^(l - e ik ^) in Eq. (3.1.13), are obtained as follows. The 
integrations over v±, t>2, v[, and v' 2 in (x^ P ) 2 produce the dipole factors (1 — e lkxn ) 

2 The exact T-matrix element (2.6.6), used to describe scattering processes at ultra-high energies 
in chapter 5, goes beyond two-gluon exchange due to the higher orders in the cosine expansion. 
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(a) (b) 




Figure 3.1: Perturbative and non-perturbative contributions to dipole-dipole scattering: 
(a) perturbative quark-quark interaction and non-perturbative (b) quark-quark, (c) string- 
string, and (d) quark-string interactions. The term quark is used genuinely for quarks and 
antiquarks. Only one diagram is shown for each type of interaction. Paired horizontal 
lines represent quark-antiquark dipoles, surrounding brackets indicate the averages over 
the dipole orientations (3.1.14), shaded areas visualize strings, curly lines denote exchanged 
perturbative gluons, and dashed and solid zig-zag lines symbolize, respectively, the non- 
perturbative non-confining and confining correlation functions. 
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and the denominator l/[(&j_^i) 2 (fc_i_r 2 ) 2 ]. This denominator multiplied with the addi- 
tional factor (r~ir~2) 2 from (x^ P ) 2 , see Eq. (3.1.5), gives the total angular dependence 




r\r\ cos 2 (0i — 02) 



{kx_r\) 2 {kx_r 2 ) 2 



{k\r\ cos 2 0i) {h\r\ cos 2 2 ) 
(cos 0i cos 02 + sin 0i sin 2 ) 



COS 2 0i COS 2 2 

j (l + 2 tan0x tan0 2 + tan 2 0! tan 2 2 ) . (3.1.17) 



± 



The first term in Eq. (3.1.17) explains the interaction between the endpoints of the 
strings with the additional factor l/kf_ in the first contribution I^ p already men- 
tioned above. The product of the second term in Eq. (3.1.17), 2tan0!tan0 2 , with 
the dipole factors vanishes after the integration over the dipole orientations 0i and 
02- The third term in Eq. (3.1.17), tan 2 x tan 2 2 , weights the different orientations 
of the strings and is characteristic for the string-string interaction (3.1.13) in our 
model. Due to this factor, the string-string interaction differs significantly from the 
interaction between the quarks and antiquarks of the dipoles known from perturba- 
tive two-gluon exchange (3.1.11). 

In Fig. 3.2 we show the perturbative integrand I p (solid line) and the non- 
perturbative integrand I NP (dashed line) of the total dipole-dipole cross section 

— * 

(3.1.10) as a function of transverse momentum \k±\ for various dipole sizes. The 
integrands have been calculated with the model parameters given in Sec. 2.7 that 
allow a good description of the experimental data of hadron and photon reactions at 
the cm. energy of ^/sq « 20 GeV with the wave functions given in Appendix B and 
the exact T-matrix element (2.6.2). Evidently, the non-perturbative integrand I NP 
governs the low momenta and the perturbative integrand I p the high momenta ex- 
changed in the dipole-dipole scattering. This behavior is a direct consequence of the 
correlation functions: The non-perturbative correlation functions (3.1.7) and (3.1.8) 
favor low momenta and suppress high momenta in comparison to the perturbative 
correlation function (2.4.4). For fixed dipole sizes \ f Dl \ and \ f D2 \ the absolute values 
of the perturbative and non-perturbative integrand are controlled, respectively, by 
the parameters (m g , M) and (G 2 , a, k) given in Sec. 2.7. With decreasing dipole 
sizes the ratio I p j I NP increases: for large dipole sizes \r Dl \ = \ f D2 \ = 1 fm, the non- 
perturbative contribution gives the main contribution to the total dipole-dipole cross 
section; for small dipole sizes \f Dl \ = \f D2 \ = 0.1 fm, the perturbative contribution 
dominates. 

As can be seen from the analytic results in Eqs. (3.1.1 1)— (3.1.16), the total 
dipole-dipole cross section (3.1.10) or the forward scattering amplitude T(so,t = 0) 
does not depend on the longitudinal quark momentum fractions z Di of the dipoles. 
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Figure 3.2: The perturbative integrand I p (solid line) and the non-perturbative integrand 
jNP gashed line) of the total dipole-dipole cross section (3.1.10) as a function of trans- 
verse momentum for various dipole sizes \f Dl \ and \ f D2 \. The integrands are shown 
for the parameters given in Sec. 2.7 that allow a good description of the experimental 
data of hadron and photon reactions at ~ 20 GeV with the wave functions given in 
Appendix B and the exact T-matrix element (2.6.2). The oscillations of the perturbative 
integrand I p originate from the Bessel function Jq(|A;_j_| \r Di \). 
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For t — the parameter z Di disappears upon the integration over Z{ since only 
\ip Di (zi, r*j)| 2 given in Eq. (3.1.2) depends on Z{. 

The structure presented for dipole-dipole scattering remains in reactions involv- 
ing hadrons and photons: the hadronic and photonic total cross sections are obtained 
from the total dipole-dipole cross section (3.1.1) by replacing \ip D .(zi, r*j)| 2 given in 
Eq. (3.1.2) with the hadron and photon wave functions given in Appendix B. As 
the total dipole-dipole cross section, the total hadronic and photonic cross sections 
are independent of the parameters which control the ^-distribution in the wave 
functions due to the normalization of the ^-distributions. The independence of the 
total hadronic cross section on the widths Azh can be seen immediately with the 
Gaussian hadron wave functions (B.l.l) 

ftMl-e^lVfc) = 1-e"^ , 

(^|tan 2 ^(l-e^)l^) = -1 + e~^ s " + J| \k ± \S h Erf 




with the error function Erf(z) = \f2pK J* dt exp(— t 2 ). These analytical results con- 
firm the Azh and z D - independence of the total dipole-proton cross section assumed 
in phenomenological models [71-73]. However, the non-forward hadronic scatter- 
ing amplitude T(so,t ^ 0) depends on the parameter Az/, as shown explicitly in 
Appendix D. Thus, the differential elastic cross section da el /dt(s,t) (5.4.1) and its 
logarithmic slope B(s,t) (5.3.1) are A^- dependent. In fact, this Azh - dependence 
is essential for the agreement with experimental data [8] as discussed in chapter 5. 

The | k± | - dependence of the perturbative and non-perturbative integrand for 
hadron-hadron, hadron-photon, and photon-photon cross sections at high photon 
virtualities is similar to the one of the perturbative and non-perturbative integrand 
of the dipole-dipole cross section shown in Fig. 3.2 for (\r Dl \ = \r D2 \ = lfm), (\f Dl \ = 
1 fm, | | = 0.1 fm), and (\r Dl \ = \r D2 \ = 0.1 fm), respectively. Of course, the 
absolute values differ and the oscillations of the perturbative integrand caused by 
the Bessel functions disappear. 
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3.2 Decomposition of the QCD String into 

Dipoles and Unintegrated Gluon Distributions 



The unintegrated gluon distribution of hadrons 3 Fhfa, k±) is a basic, universal quan- 
tity convenient for the computation of many scattering observables at small x. It 
is the central object in the BFKL [40] and CCFM [41] evolution equations. Upon 
integration over the transverse gluon momentum \k±\ it leads to the conventional 
gluon distribution xGh{x,Q 2 ) used in the DGLAP evolution equation [42]. The 
unintegrated gluon distribution is crucial to describe processes in which transverse 
momenta are explicitly exposed such as dijet [38] or vector meson [39] production 
at HERA. Its explicit \k±\ dependence is particularly suited to study the interplay 
between soft and hard physics. In this section an exact representation of the string 
as a collection of stringless quark-antiquark dipoles is presented that allows us to 
extract the perturbative and non-perturbative contributions to !Fh( x , k±) from our 
total dipole-hadron cross section via \k±\ -factorization. 

We calculate the unintegrated gluon distribution !Fh(x, k 2 ) with the T-matrix 
element in the limit of small X"f unc tions (2.6.3). Following Sec. 2.6, we give a 
strong energy dependence to the perturbative contribution \ F and a weak one to 
the non-perturbative contribution x NF \ see also Eqs. (2.6.4) and (2.6.5), 

NP 

(x Np r - (x Np r(^) (3.2.i) 

where the values of the exponents e p and e NP are given in Sec. 2.7 and xq = 2.4 x 10 -3 
is adjusted to reproduce at Q 2 = 1 GeV 2 the integrated gluon distribution of the 
proton xG p (x, Q 2 ) extracted from the HERA data [79]. The small- x limit (2.6.3) 
considered here is applicable only for x = Q 2 /s > 10~ 4 . For x < 10~ 4 , the full 
T-matrix element (2.6.2) has to be used which satisfies the S'-matrix unitarity and 
leads to a successful description of many reactions as shown in chapter 4 and 5. The 
structure of the interactions worked out in this section, however, is independent of 
the phenomenological energy dependence considered. 

The total dipole-hadron cross section o~E>h(x, \ f D \) is obtained from the total 
dipole-dipole cross section (3.1.1) by replacing \i^ D2 {z 2) r 2 ) \ 2 with a squared hadron 
wave function 1^,(^2, ^2) | 2 - Accordingly, the x-dependent total dipole-hadron cross 

3 The word hadron and the subscript h is used genuinely for hadrons and photons in this section. 
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section reads 



ODhfa \r D \) = J dk\ 



(3.2.2) 



x 



(4na s (k 2 ± )Y{ i&^ikl) ' ( 1 - J (\k ± \ \r D \) ) (^|1 - e ik ^ |^> \ ( ^ 



x 



7T^ 2 

24 



1 2 



l-M\k±\\r D \)) Ml - e ik ^\ij> h ) 



+ 



-1+ ^(-1; 1 1; ) ty fc | tan 2 2 (l - e^)|^> 



x / 



with the Bessel function Jo(|^i||^d|) an d the generalized hypergeometric function 
iF 2 (— 1/2; 1/2, 1; — kfr^/A) derived in the previous section. 

For dipole sizes \r D \ — > 0, the perturbative contribution to cr Dh (x, \r D \) is known 
to vanish quadratically with decreasing dipole size, o~Dh( x , \r D \) oc r 2 . This behavior 
reflects the weak absorption of a small color-singlet dipole in the hadron and is 
known as color transparency. It can be seen immediately from Eq. (3.2.2) as 



i-M\k±\K\) 



± D for I fx, I — > and finite \k± 



(3.2.3) 



The non-perturbative contribution to cr D/l (x, \ f D \) gives color transparency as well 
since the generalized hypergeometric function behaves as 



1 1 —k 2 r 2 
-l + 1 F 2 (--;-,l;^^ 



k\rl 



for |fx,| — > and finite \k±\. (3.2.4) 



For large dipole sizes, \r D \ > 1 fm, the perturbative contribution to aoh( x , 1^1) 
describing interactions of the quark and antiquark of the dipole with the hadron 
saturates since 

(l - M\k±\\fb\)) « 1 for large |fc ± ||fb| • (3.2.5) 

In contrast, the non-perturbative contribution to a£>h{x, \f D \) increases linearly with 
increasing dipole size, ach( x , \f D \) oc \r D \. This linear increase is generated by the 
interaction of the string of the dipole with the hadron: The string elongates linearly 
with the dipole size \ f D \ and, thus, has a linearly increasing geometric cross section 
with the hadron. Indeed, this feature of the string can be seen analytically since 

1 + iF 2 {~; \, 1; Z ^)) oc \k ± \\r D \ for large |4||r D | ■ (3.2.6) 
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When considered in Euclidean space-time, the same string gives also the linear con- 
fining potential between a static quark and antiquark at large qq separations [17, 31] 
as shown briefly at the end of Sec. 2.4. Thus, the behavior of the total dipole-hadron 
cross section is related to the confining potential. Furthermore, as we are working in 
the quenched approximation, there is no string breaking through dynamical quark- 
antiquark production at large dipole sizes. String breaking is expected to stop the 
linear increase of the total dipole-hadron cross section at dipole sizes of \ f D \ > 1 fm 
analogous to the saturation of the static qq potential seen for large qq separations 
on the lattice in full QCD [26, 80]. 

In Fig. 3.3 we show the perturbative (solid line) and non-perturbative (dashed 
line) contributions to the total dipole-proton cross section (Td p {x, \r D \) as a function 
of the dipole size \r D \ at x = Xq = 2.4 • 10~ 3 , where the perturbative contribution is 
multiplied by a factor of 10. The dipole-proton cross section is computed with the 
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Figure 3.3: Perturbative (solid line) and non-perturbative (dashed line) contributions to 
the total dipole-proton cross section an P (x, \f D \) as a function of the dipole size \f D \ at 
x = xq = 2.4 • 10 -3 . with the perturbative contribution multiplied by a factor of 10. The 
perturbative contribution shows color transparency at small \r D \ and saturates at large 
\f D \. The non-perturbative contribution shows also color transparency at small \f D \ but 
increases linearly with increasing \r D \. 
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simple Gaussian proton wave function (B.l.l) and illustrates the general features 
discussed above: The perturbative contribution shows color transparency at small 
\r D \ and saturates at large \f D \. The non-perturbative contribution shows also color 
transparency at small \ f D \ but increases linearly with increasing \f D \. 

Our result (3.2.2) shows that the \k±\- dependence of the hadron constituents fac- 
torizes from the rest of the process in both the perturbative and the non-perturbative 
contribution to the total dipole-hadron cross section. This factorization - known 
in perturbative QCD as \k±\ -factorization [81] - allows to define the unintegrated 
gluon distribution J r / l ( follows [69, 72, 82] 

4 2 2 r [1 - Jo{\k±\\r D \)) 
ob h (x,\? D \) = — ± W V — — - 1 a s {kl)T h {x,kl) . (3.2.7) 
6 J {\k±\\r D \y 

For small dipole sizes \ f D \, this equation together with (3.2.3) and the integrated 
gluon distribution 

rQ 2 

xG h (x,Q 2 ) = / dklT h (x,kl) (3.2.8) 
Jo 

leads to the widely used perturbative QCD relation [83] 

2 2 

a Dh (x, \r D \) = ^ [a s (Q 2 )xG h (x,Q 2 )} Q2 = c/r , , (3.2.9) 

where c ~ 10 is estimated from the properties of the Bessel function Jod^il I'd I) [84]. 

To extract the unintegrated gluon distribution J-h(x, kf), we compare (3.2.7) 
with (3.2.2) using the following mathematical identity 

(-l + iF 2 (-l;l 1;^)) = J\ ^ {l - Jo(\k ± mO) ■ (3-2.10) 

As discussed in the previous section, the lhs of (3.2.10) results from the string aver- 
aged over all orientations (3.1.16). Thus, the string confining the quark- ant iquark 
dipole of length |f^| can be represented as an integral over stringless dipoles of sizes 
£\r D \ with < £ < 1 and a dipole number density of n(£) = l/£ 2 . As visualized in 
Fig. 3.4, the string-hadron scattering process reduces to an incoherent superposition 
of stringless dipole-hadron scattering processes with dipole sizes < £\r D \ < \f D \ 
and dipole number density n(£) = l/£ 2 . The decomposition of the string into many 
smaller stringless dipoles via (3.2.10) behaves similar to the wave function of a qq 
onium state in the large N c limit [36,37]: The numerous gluons emitted inside the 
onium state can be considered as many qq dipoles in the large N c limit. 
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Figure 3.4: The string of length \r D \ is made up of stringless dipoles of size ^\r D \ with 
< £ < 1 and dipole number density n(£) = l/£ 2 . The string-hadron scattering process 
reduces to an incoherent superposition of stringless dipole-hadron scattering processes. 



Inserting (3.2.10) into (3.2.2) and rescaling the momentum variable \k[_\ = £\k±\. 
the string-hadron (sh) contribution to the total dipole-hadron cross section (3.2.2) 
becomes 



<Mrn\) = ~ j dk>* (l - U\ki L \\r D \) 



(3.2.11) 
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24 



2 /•! 
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(^|tan 2 ^(l-e^/^)|^> 



x / 



The dipole factor (1 — J (IM_| |^d|) indicates that the string-hadron interaction has 
been rewritten into a stringless dipole-hadron interaction. The string confining the 
dipole has been shifted into the hadron. Comparing (3.2.7) with our new expression 
for the total dipole-hadron cross section, 



<TDh(x, \r D \) 



9 4vr 



dkl [l-M\k L \\r D \) 
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(^|tan 2 </) 2 (l-e^/^)|^) 
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(3.2.12) 
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one obtains the unintegrated gluon distribution 

k 2 



6n 3 a s (kD 



x 



+ 



„-n'(2)n.2 
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7T G 2 
24 

2 fl 







'1-/0 



i£>i (2) (lfc^ 



X 
2 



(^|tan^ 2 (l-e^/^)|^) 



A |l-e^|^ fc > 

^Y" . (3.2.13) 

x J 



This result shows explicitly the microscopic structure of the perturbative and non- 
perturbative contribution to J-"h(x, k 2 ). It is valid for any hadron wave function. To 
get numerical values for the unintegrated gluon distribution !Fh(x, kf_), the hadron 
wave functions must be specified. 



3.3 Numerical Results for Unintegrated Gluon 
Distributions in Hadrons and Photons 

In this section we present the unintegrated gluon distribution (3.2.13) for protons, 
pions, kaons, and photons computed with the Gaussian hadron wave function (B.l.l) 
and the perturbatively derived photon wave functions (B.2.1) and (B.2.2). To ac- 
count for non-perturbative effects at low photon virtuality Q 2 in the photon wave 
functions, quark masses rrif(Q 2 ) are used that interpolate between the current quarks 
at large Q 2 and the constituent quarks at small Q 2 [20] as discussed in Appendix B. 

The unintegrated gluon distribution of the proton !F p (x, kf) is shown as a function 
of transverse momentum \kj_\ at x = 10" 1 , 10~ 2 , 10~ 3 , and 10~ 4 in Fig. 3.5. Fig- 
ure 3.6 illustrates the interplay of the perturbative (solid line) and non-perturbative 
(dashed line) contributions to |A;j_|^( function of transverse momentum 

\kj_\ for the same values of x. At small momenta \kj_\, the unintegrated gluon dis- 
tribution is dominated by the non-perturbative contribution and behaves as l/|fcj_|. 
This behavior is a string manifestation as it reflects the linear increase of the to- 
tal dipole-proton cross section at large dipole sizes. In contrast, the saturation 
model of Golec-Biernat and Wiisthoff [72] shows the behavior J-^ BW (x, kf) oc k]_ for 

small momenta. With increasing \k±\, the non-perturbative contribution to !F p (x, kf_) 
decreases rapidly which results from the strong suppression of large momenta by 
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the non-perturbative correlation functions iD'i \k±) and iD^ 2 \kf] given in (3.1.7) 
and (3.1.8). For \k±\ > 1 GeV, the perturbative contribution dominates the unin- 
tegrated gluon distribution. It drops as 1/k 2 in accordance with the perturbative 
correlation function it)'p{kf) given in (3.1.6). This perturbative QCD result is 
not reproduced by the phenomenological model of Golec-Biernat and Wiisthoff [72] 
which predicts a Gaussian decrease of T v {x, kj_) with increasing \k±\. 




10 -4 I I I I I I I 

1 2 3 4 5 6 
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Figure 3.5: The unintegrated gluon distribution of the proton ^,(x,/cj) as a function of 
transverse momentum at Bjorken-x values of 10 , 10 -2 , 10 -3 and 10 -4 . 

The x- dependence of J-p(x, kf_) can be seen in Figs. 3.5 and 3.6: With decreasing 
x, the perturbative contribution increases much stronger than the non-perturbative 
contribution which results from the energy exponents e p ^> e NP in (3.2.1) neces- 
sary to describe the experimental data within our model as discussed in Sec. 2.6. 
Moreover, the perturbative contribution extends into the small- \k±\ region as x 
decreases. Indeed, the soft-hard transition point moves towards smaller momenta 
with decreasing x as shown in Fig. 3.6. Such a hard-to-soft diffusion is observed also 
in [85] where the unintegrated gluon distribution has been parametrized to repro- 
duce the experimental data for the proton structure function F 2 (x, Q 2 ) at small x. 
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Figure 3.6: The unintegrated gluon distribution of the proton F p {x, k±) times the tran- 
verse momentum \k±\ as a function of \k±\ at Bjorken-x values of 10 _1 , 1CT 2 , 10~ 3 and 
10" 4 . 



The opposite behavior is obtained in the approach of the color glass condensate [86]: 
With decreasing x, gluons are produced predominantely in the high - \kj_\ region of 
lower density and weaker repulsive interactions. 

In Fig. 3.7, the unintegrated gluon distribution of the proton, pion, and kaon 
Thix, kf) times the transverse momentum \k±\ is shown as a function of \k±\ at 
x = 10 -3 . The hadrons are characterized by different values for Azh and Sh in 
the hadron wave function (B.l.l). However, ^(x, /c 2 ) depends only on Sh- Due to 
the normalization of the hadron wave functions, Azh disappears upon the integra- 
tion over Z{ as can be seen directly from (3.1.18) and (3.1.19). At small momenta, 
J r h(x,k^) oc Si is found from (3.1.18), (3.1.19), and (3.2.13). It becomes visible 
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\k ± \ [GeV] 

Figure 3.7: The unintegrated gluon distribution of the proton, pion, and kaon J-h(x,k±) 
times the transverse momentum \k±\ as a function of \k±\ at Bj or ken- variable x = 10~ 3 . 

in Fig. 3.7 for the chosen hadron extensions: S p = 0.86 fm, S w = 0.607 fm, and 
Sk = 0.55 fm. At large momenta where the perturbative contribution dominates, 
the dependence on Sh vanishes as can be seen from (3.1.18) and the unintegrated 
gluon distributions of protons, pions, and kaons become identical. Of course, this 
behavior results from the wave function normalization being identical for protons, 
pions, and kaons with two valence constituents which are the quark and antiquark 
in the pion and kaon and the quark and diquark in the proton. At large \k±_\, i.e., 
high resolution, the realistic description of protons as three-quark systems becomes 
necessary. In fact, the three-quark description of protons leads to a different per- 
turbative contribution in (3.2.13): the quark counting factors of 2 (appropriate for 
mesons) in the square brackets in (3.1.11) are substituted by factors of 3 (appropri- 
ate for baryons) [34,35]. At other values of x, the unintegrated gluon distributions 
of protons, pions and kaons show the same features. 

Photons are particularly interesting because the transverse size of the quark- 
antiquark dipole into which a photon fluctuates is controlled by the photon virtual- 
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Figure 3.8: The unintegrated gluon distribution of the transverse polarized photon 
J-^(x,k^_,Q 2 ) times the tranverse momentum as a function of \k±\ at photon vir- 
tualities of Q 2 = 0, 1, 10, and 100 GeV 2 and Bj or ken- variable x = 10 -3 . 

ity Q 2 (cf. Appendix B) 

2 

where m u (Q 2 ) is the running w-quark mass given in Appendix B. In Figs. 3.8 
and 3.9, the unintegrated gluon distribution of transverse (T) and longitudinal (L) 
photons L (x, kj_, Q 2 ) times the transverse momentum \k±\ is shown as a function 

of \k±\ for various photon virtualities Q 2 at x = 10~ 3 . With increasing Q 2 , i.e., 
decreasing "photon size" |r 7 * |, the ratio of the perturbative to the non-perturbative 
contribution to J^,* (x, k 2 , Q 2 ) increases. This behavior has already been discussed 
in Sec. 3.1 on the level of dipole-dipole scattering, where decreasing dipole sizes 
increase the ratio of the perturbative to non-perturbative contribution, see Fig. 3.1. 
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Figure 3.9: The unintegrated gluon distribution of the longitudinally polarized photon 
T^* L (x, k 2 , Q 2 ) times the tranverse momentum \k±\ as a function of \k±\ at photon virtual- 
ities of Q 2 = 0.01, 1, 10, and 100 GeV 2 and Bj or ken- variable x = 10" 3 . 

Due to the different Q 2 - dependence in the transverse and longitudinally polarized 
photon wave functions, F^{x, k]_, Q 2 ) decreases continously with increasing Q 2 while 
J~il{ x -i k±, Q 2 ) increases for Q 2 < 1 GeV 2 and decreases for Q 2 > 1 GeV 2 . 

In Fig. 3.10, the unintegrated gluon distributions of hadrons and photons dis- 
cussed above are compared. The unintegrated gluon distribution of real (Q 2 = 0) 
photons (|r 7T | ~ Sh) is suppressed by a factor of order a in comparison to the one 
of the hadrons otherwise its shape is very similar. The suppression factor comes 
from the photon- dipole transition described by the photon wave functions given 
in Appendix B. The ratio of the perturbative to the non-perturbative contribution 
of unintegrated gluon distributions increases as one goes from hadrons to photons 
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Figure 3.10: The unintegrated gluon distribution of the proton, pion, kaon, and transverse 
and longitudinally polarized photon J-h,-y* (x, k 2 ) times the transverse momentum \ k±\ as 

a function of at Bj or ken- variable x = 10 -3 . Results for transverse polarized photons 
are shown for photon virtualities of Q 2 = and 10 GeV 2 and results for longitudinally 
polarized photons for photon virtualities of Q 2 = 0.01 and 10 GeV 2 . 

with high virtuality Q 2 . Since the wave functions of protons, pions, and kaons are 
normalized to the same value, the unintegrated gluon distributions of these hadrons 
are, as mentioned above, identical at large \k±\ and do not depend on the hadron 
size. In contrast, the Q 2 - dependence of the photon wave functions leads to Q 2 - 
dependent, i.e., "photon size" - dependent, unintegrated gluon distributions at large 
With increasing \k±\, the unintegrated gluon distributions of hadrons and pho- 
tons become parallel in line with the vanishing dependence on the specific form of 
the wave function. 

In Fig. 3.11, the integrated gluon distribution of the proton xG p (x, Q 2 ) is shown 
as a function of Bjorken-x at photon virtualities of Q 2 = 1, 5, and 20 GeV 2 . Recall 
that the parameter x = 2.4 • 10~ 3 has been adjusted in the previous section such 
that the experimental data of xG p (x, Q 2 ) at Q 2 = 1 GeV 2 [79] are reproduced. For 
x > 10~ 3 , xG p (x, Q 2 ) is mainly determined by non-perturbative physics as can be 
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Figure 3.11: The gluon distribution of the proton xG p (x, Q 2 ) as a function of Bjorken-x 
at photon virtualities of Q 2 = 1, 5, and 20GeV 2 . 

seen from Figs. 3.5 and 3.6. Perturbative physics becomes relevant for x < 10~ 3 
and generates the steep increase of xG p (x, Q 2 ) with decreasing x at fixed Q 2 . Also 
the rise of xG p (x, Q 2 ) with increasing Q 2 at fixed x results from the perturbative 
contribution. For x <C 10~ 4 , we show explicitly in Sec. 4.4 that multiple gluonic 
exchanges contained in the full T-matrix element (2.6.6) slow down the powerlike 
increase of xG p (x, Q 2 ) with decreasing x in accordance with S'-matrix unitarity 
constraints. 
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3.4 Comparison with Other Work 



In this section, we compare the unintegrated gluon distribution of the proton ex- 
tracted from our loop-loop correlation model (LLCM) with those obtained from the 
saturation model of Golec-Biernat and Wusthoff (GBW) [72], the derivative of the 
Gliick, Reya, and Vogt (GRV) parametrization of xG p (x, Q 2 ) [87], and the approach 
of Ivanov and Nikolaev (IN) [85]. 



In the approach of Golec-Biernat and Wusthoff [72], the unintegrated gluon 
distribution of the proton is extracted from the total dipole-proton cross section by 
inverting (3.2.7), 

F p (x,kl) = -^Rl(x)klexp{-R 2 (x)kl) with R 2 Q (x) = -L (fV (3 .4 .1) 

4 71 QL S LJq \Xq J 

where the parameter a = 29.12 mb, a s = 0.2, Q = 1 GeV, A = 0.277, and x = 
0.41-10 -4 are obtained from a fit to the proton structure function F 2 (x, Q 2 ) including 
charm quarks in the photon wave function. Note, however that the GBW approach 
uses the dipole-proton cross section of the GBW model on the lhs of (3.2.7) which 
implies multiple gluon exchanges while the rhs of (3.2.7) describes only two-gluon 
exchange as discussed in Sec. 3.1. Moreover, as demonstrated below, the large kj_- 
behavior of the unintegrated gluon distribution (3.4.1) deviates significantly from the 
DGLAP results. This mismatch motivated the recent modifications of the saturation 
model [88]. 



The unintegrated gluon distribution of the proton is also computed from the 
integrated gluon distribution xG p (x, Q 2 ) by inverting (3.2.8) 



J- p (x, fcj_) 



dxG p (x, Q 2 



dQ 2 



(3.4.2) 



Here, we use for the integrated gluon distribution xG p (x, Q 2 ) the leading order (LO) 
parametrization of Gliick, Reya, and Vogt (GRV) [87], which covers the kinematic 
region 10" 9 < x < 1 and 0.8 GeV 2 < Q 2 < 10 6 GeV 2 . 



Ivanov and Nikolaev [85] have constructed a two-component (soft + hard) ansatz 
for the unintegrated gluon distribution of the proton 

k 2 2 kl 

J- {X, kj_) = J 7 30 ft kj_) 2 + ^liard(^) kj_) , 2 , 7,2 (3.4.3) 
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with the soft and hard component 



^ r so ft(fci) 

J~ha,rd{x, kj_~) 
-(B) 



Osoft Cf N c 



a s {kl 
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3kl 
A 2 



(3.4.4) 



where .7-"^ (fc 2 ) has the same form as (3.4.4) with the parameters a pt and /i pt in- 
stead of a so ft and /i so f t and T p t(x,kf) is the derivative of the integrated gluon dis- 
tribution (3.4.2). In the IN approach, the running coupling is given by a s (kf) = 
min{0.82, (4tt)/(/3 logf^/A^])}. With the GRV-parametrization [87] for the 
integrated gluon distribution, the structure function of the proton F 2 (x,Q 2 ) has 
been described successfully using the following parameters: k 2 s 



(1 + 0.0018 log(l/x) 4 ) - 5 
Ql = 0.895 GeV 2 , # 



. Osoft — 2, 

9, and A QC d 



l pt — 1) /^soft 

= 0.2 GeV. 



0.1 GeV, 



3 GeV 2 , k\ = 
0.75 GeV, 



In Fig. 3.12, we show the LLCM, GBW, GRV, and IN results for the unintegrated 
gluon distribution of the proton !Fp(x, kf) as a function of transverse momentum 
squared k\ for x = 10 _1 , 10 -2 , 10 -3 , and 10~ 4 . At small transverse momenta, 
kl < 0.1 GeV 2 , our model gives the largest values for J- p (x, kf_). As mentioned in 
the previous section, our LLCM unintegrated gluon distribution increases as 1/ a/AJ 
with decreasing k\ as a consequence of the linear increase of the total dipole-proton 
cross section at large dipole sizes. In contrast, for k\ — > 0, the unintegrated gluon 
distribution of GBW decreases as k\ and the one of IN as kf_. In the perturbative 
region, k\ > 1 GeV 2 , the unintegrated gluon distribution of the LLCM becomes 
smaller than the one of GRV and IN but is still larger than the one of GBW. 
Moreover, the LLCM, GRV, and IN unintegrated gluon distributions become parallel 
for x < 10~ 2 and drop as 1/A; 2 for large k\. This perturbative QCD behavior 
is not reproduced by the GBW unintegrated gluon distribution which decreases 
exponentially with increasing k\. 



The x-dependence of the LLCM, GBW, GRV, and IN unintegrated gluon dis- 
tributions J- p (x, kf) is shown for transverse momenta squared k\ = 0.1, 0.5, 10, 
and 50 GeV 2 in Fig. 3.13. In the non-perturbative region, kl = 0.1 GeV 2 and 
k]_ = 0.5 GeV 2 , the LLCM, GBW, and IN unintegrated gluon distributions show a 
weak x- dependence. In the perturbative region, k 2 = 10 GeV 2 and k\_ = 50 GeV 2 , 
the x-dependence of the unintegrated gluon distributions becomes stronger. The 
LLCM, GRV, and IN unintegrated gluon distributions show nearly the same rise 
with decreasing x. In contrast, the GBW unintegrated gluon distribution increases 
much faster decreases. 
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Figure 3.12: The unintegrated gluon distribution of the proton T p (x,k]_) as a func- 
tion of transverse momentum squared kj_ at Bjorken-x values of 10 _1 , 10~ 2 , 10~ 3 , and 
10~ 4 . The different curves are obtained from our loop-loop correlation model (LLCM), 
the Golec-Biernat and Wiisthoff (GBW) model [72], the derivative of the Gliick, Reya, 
and Vogt (GRV) parametrization of xG p (x,Q 2 ) [87], and the Ivanov and Nikolaev (IN) 
approach [85]. 
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Figure 3.13: The unintegrated gluon distribution of the proton F p (x, kj) as a function of 
Bjorken-x at transverse momenta squared k]_ = 0.1, 0.5, 10, and 50GeV 2 . The different 
curves are obtained from our loop- loop correlation model (LLCM), the Golec-Biernat 
and Wiisthoff (GBW) model [72], the derivative of the Gluck, Reya, and Vogt (GRV) 
parametrization of xG p (x, Q 2 ) [87], and the Ivanov and Nikolaev (IN) approach [85]. Note 
that the GRV parametrization is only available for k]_ > 0.8 GeV 2 . Moreover, the GBW 
result for ^(x.fcf) is below 10" 5 for k\ = 50 GeV 2 . 
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In addition to the one-scale unintegrated gluon distributions, T(x, kj), discussed 
in this work, there exist also two-scale unintegrated gluon distributions, T{x, k]_, /i 2 ). 
In the CCFM evolution equation [41], this additional scale /i 2 is related to the max- 
imum angle allowed in the gluon emission. Two-scale unintegrated gluon distribu- 
tions are obtained in the approach of Blumlein [89], Jung and Salam [90], Kimber, 
Martin, and Ryskin [91], and in the linked dipole chain (LDC) model [92,93]. A 
comparison of their results can be found in [93,94] where also the one-scale unin- 
tegrated gluon distributions of Kwiecinski, Martin and Stasto [95] , and Ryskin and 
Shabelski [96] are discussed. 
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Chapter 4 



Impact Parameter Profiles 
and Gluon Saturation 



In this chapter we study saturation effects that manifest the unitarity of the S- 
matrix at high cm. energies. We investigate the scattering amplitude in impact 
parameter space where the S-matrix unitarity limit becomes most explicit. In fact, 
the S'-matrix unitarity imposes the black disc limit on the height of such impact 
parameter profiles. We show explicitly that our model respects the black disc limit. 
Furthermore, the width of the impact parameter profiles is shown to increase loga- 
rithmically at asymptotic energies as needed to guarantee the Froissart bound [7]. 

We compute the impact parameter profiles for hadron-hadron and longitudi- 
nal photon-proton scattering. Concrete energy values are determined at which the 
profiles saturate at the black disc limit for small impact parameters. The impact 
parameter profiles provide an intuitive geometrical picture for the energy depen- 
dence of the scattering process as they illustrate the evolution of the opacity and 
size of the interacting particles with increasing cm. energy. We use the explicit 
black disc limit of the profiles in the following chapter to localize saturation effects 
in experimental observables. 

Using a leading twist, next-to-leading order DGLAP relation, we estimate the 
impact parameter dependent gluon distribution of the proton xG(x,Q 2 , \b±\) from 
the profile function for longitudinal photon-proton scattering. We find a low-x 
saturation of xG(x,Q 2 , \b±\) as a manifestation of the S'-matrix unitarity. The im- 
plications on the integrated gluon distribution xG(x, Q 2 ) are studied and compared 
with complementary investigations of gluon saturation. 

The results shown in this chapter are obtained with the T-matrix element given 
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in (2.6.6), the model parameters explained in Sec. 2.7, and the hadron and photon 
wave functions discussed extensively in Appendix B. 



4.1 S- Matrix Unitarity and Impact Parameter Profiles 

In this section we show that our model respects the S'-matrix unitarity limit in 
impact parameter space of the scattering amplitude. 

The impact parameter dependence of the scattering amplitude is given by T(s, \b±_\), 

T(s,t = -tf) = 4:sjd 2 b ± e liLk T(s,\b ± \) (4.1.1) 

and in particular by the profile function 

J(s,\b x \) = 2ImT(s,\b ± \) , (4.1.2) 

which describes the blackness or opacity of the interacting particles as a function 
of the impact parameter \b±\ and the cm. energy yfs. In fact, the profile func- 
tion (4.1.2) determines every observable if the T- matrix is — as in our model - 
purely imaginary. 

The S'-matrix unitarity, SS^ = S^S = 11, leads directly to the unitarity condition 
in impact parameter space [97, 98] 

ImT(s, |6l|) = \T(s, \b ± \)\ 2 + G inel (s, \b ± \) , (4.1.3) 

where G ine i(s, \b±\) > is the inelastic overlap function [99]. 1 This unitarity condi- 
tion imposes an absolute limit on the profile function 

< 2\T(s,\b ± \)\ 2 < J(s,\b ± \) < 2 (4.1.4) 

and the inelastic overlap function, G ine i(s, \b±\) < 1/4. At high energies, however, 
the elastic amplitude is expected to be purely imaginary. Consequently, the solution 
of (4.1.3) reads 

J(s, \b ± \) = 1 ± ^1-4^(3,1^1) (4.1.5) 

and leads with the minus sign corresponding to the physical situation to the reduced 
unitarity bound 

< J(s, \b x \) < 1 . (4.1.6) 

integrating (4. 1 .3) over the impact parameter space and multiplying by a factor of 4 one obtains 
the relation a tot {s) = a el {s) + a mel {s). 
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Reaching the black disc limit or maximum opacity at a certain impact parameter 
\b±\, J(s, \b±\) = 1, corresponds to maximal inelastic absorption Gi ne i(s, \b±\) = 1/4 
and equal elastic and inelastic contributions to the total cross section at that impact 
parameter. 

In our model, every reaction is reduced to dipole-dipole scattering with well 
defined dipole sizes |r*j| and longitudinal quark momentum fractions Z{. The unitarity 
condition in our model becomes, therefore, most explicit in the profile function 



Jdd(s, \b ± \,zi, |fi|,z 2 , \f 2 \ 



2tt 



2vr L 



S D D(s,b ± ,z 1 ,r 1 ,z 2 ,r< 2 ) , (4.1.7) 



where 4>i describes the dipole orientation, i.e. the angle between fi and b±, and Sdd 
describes elastic dipole-dipole scattering 



S 



DD 



COS 



cos 



cos 



';X NP (s, 



cos 



X P (s)) (4.1.8) 



with the purely real-valued eikonal functions x ( s ) an d X ( s ) defined in (2.6.4). 
Because of \Sdd\ < 1, a consequence of the cosine functions in (4.1.8) describing 
multiple gluonic interactions, Jdd respects the absolute limit (4.1.4). Thus, the 
elastic dipole-dipole scattering respects the unitarity condition (4.1.3). At high 
energies, the arguments of the cosine functions in Sdd become so large that these 
cosines average to zero in the integration over the dipole orientations. This leads 
to the black disc limit J D % D X = 1 reached at high energies first for small impact 
parameters. 



If one considers the scattering of two dipoles with fixed orientation, the inelastic 
overlap function obtained from the unitarity constraint (4.1.3), 



C DD i 



1 - 



cos 



:X NP (s] 



cos 



;X P (s: 



cos 




shows nonphysical behavior with increasing energy. This behavior is a consequence 
of aritifically fixing the orientations of the dipoles. If one averages over the dipole 
orientations as in all high-energy reactions considered in this work, no unphysical 
behavior is observed. 



4.2 Profile Function for Hadron-Hadron Scattering 



In this section we show the cm. energy and impact parameter dependence of the 
hadron-hadron scattering amplitude, determine the energy values at which satura- 
tion effects set in as a manifestation of the S'-matrix unitarity, and compute the 
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logarithmic rise of the profile width at asymptotic cm. energies which is equivalent 
to the Froissart bound [7]. 



The profile function for hadron-hadron (hh) scattering 
Jhh(s,\b±\)= \dz x d 2 r 1 dz 2 d 2 r 2 \ip h (z h n)\ 2 \ip h (z2, r 2 )\ 2 1 - Sdd(s, b±, zi, n, z 2 , r 2 ) 



(4.2.1) 

is obtained from (4.1.7) by weighting the dipole sizes \ fi\ and longitudinal quark mo- 
mentum fractions Zi with the hadron wave functions \ij)h(zi, fi)\ 2 from Appendix B. 

The profile functions for proton-proton, pion-proton, and kaon-proton scattering 
are shown in Fig. 4.1 for cm. energies y/s from 10 GeV to 10 8 GeV. The pro- 
files are obtained with model parameters from Sec. 2.7 and the hadron extensions 
S p = 0.86 fm, S n = 0.607 fm, and Sk = 0.55 fm used in the Gaussian hadron wave 
functions. Up to yfs ~ 100 GeV, the profiles have approximately a Gaussian shape. 
Above y^i = 1 TeV, they significantly develop into broader and higher profiles until 
the black disc limit is reached for \fs ~ 10 6 GeV and = 0. At this point, the 
cosine functions in Sdd average to zero 

jd Zl d 2 n jdz 2 d 2 r 2 \ij h {z l ,r l )\ 2 \ij h {z 2 ,r 2 )\ 2 S DD {^~s> 10 6 GeV, 1^1 =0,...) ^0 

(4.2.2) 

so that the hadron wave function normalizations determine the maximum opacity 



J™ ax = dz l d 2 r l 



iJdz 2 d 2 r 2 \Mzi,n)\ 2 \M^r2)\ 2 = l • (4.2.3) 



Once the maximum opacity is reached at a certain impact parameter, the profile 
function saturates at that and extends towards larger impact parameters with 
increasing energy. Thus, the multiple gluonic interactions important to respect the 
S"-matrix unitarity constraint (4.1.3) lead to saturation for > 10 6 GeV. The 
above behavior of the profile functions illustrates the evolution of the hadron with 
increasing cm. energy. The hadron is gray and of small transverse size at small 
\/s but becomes blacker and more transversally extended with increasing y 7 ? until 
it reaches the black disc limit in its center at ^/s ~ 10 6 GeV. Beyond this energy, 
the hadron cannot become blacker in its central region but in its periphery with 
continuing transverse growth. Furthermore, the hadron boundary seems to stay 
diffusive as claimed also in [100]. 

The differences between the profile functions for pp (pp) , 7r ± p, and K^p scattering 
especially at low energies result from the different transverse hadron extensions, 
S p = 0.86 fm > S n = 0.607 fm > Sk = 0.55 fm, cf. Appendix B. For the smaller 
pion and kaon extensions the short distance physics described by the perturbative 
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Figure 4.1: The profile function for proton-proton, pion-proton, and kaon-proton scat- 
tering Jhh(s,\b±\) is shown versus the impact parameter for cm. energies y/s from 
10 GeV to 10 8 GeV. The unitarity limit (4.1.4) corresponds to Jhh(s,\b±\) = 2 and the 
black disc limit (4.1.6) to Jhh{s-> \b±\) = 1- 
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component becomes more important and leads to a stronger energy increase (e p = 
0.73 > e = 0.125) of the % ± p and K ± p profile functions as compared to the one 
of the pp (pp) scattering in the cm. energy range 10 GeV < ^/s < 10 6 GeV. The 
black disc limit (4.2.3) is identical for up, Kp, and pp (pp) scattering due to the 
same pion, kaon, and proton wave function normalization (B.1.2). The smaller size 
of pions and kaons in comparison to the one of protons demands, however, slightly 
higher cm. energies to reach the black disc limit and yields narrower n ± p and K^p 
profile widths for energies somewhat larger than a/s = 10 6 GeV. At asymptotic 
energies, however, the profile widths become universal (hadron size independent) as 
shown below. 

According to our model the black disc limit will not be reached at LHC. Our 
prediction of a/s « 10 6 GeV = 10 3 TeV for the onset of the black disc limit in 
proton-proton collisions is about two orders of magnitude beyond the LHC energy 
a/s = 14 TeV. This is in contrast, for example, with [101], where the value predicted 
for the onset of the black disc limit is \/s = 2 TeV, i.e. small enough to be reached at 
LHC. However, we feel confidence in our LHC prediction since our profile function 
Jpp(s, \b±\) yields good agreement with experimental data for cross sections up to 
the highest energies as shown in chapter 5. 

The discussion in the previous section, the black disc limit in Eq. (4.2.3), and 
Fig. 4.1 show that our model respects the iS-matrix unitarity in impact parameter 
space of the scattering amplitude. This, however, is not enough to guarantee the 
Froissart bound [7] for total hadronic cross sections at asymptotic cm. energies. 
Since the latter are obtained by integrating the impact parameter profiles over all 
impact parameters, difficulties with the Froissart bound can arise if the black disc 
radius R(s), i.e., the impact parameter range for which Jhh(s, \b±\) ~ 1, increases 
with growing cm. energy faster than ln(s). A calculation of the energy dependence 
of R(s), which agrees with the numerical results shown in chapter 5, shows below 
that the Froissart bound is respected in our model. 

For impact parameters much larger than the dipole sizes involved in the interac- 
tion, \b±\ ^> |ri|, |r 2 |, the leading term of the x P -function (2.5.18) which is needed 
in the computation of R(s) reads 



X P (s,fi,r 2 A) 




m G \b±\ - 



m,Q \bj_ 



\4|M 2 2\b±\J 




5 



2m G 



+ m G \fl\\f 2 \ cos(0i) cos(0 2 ) 



+ 



16. 



416, I 2 
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(4.2.4) 



where denotes the angle betweeen the dipoles fi and r*2 and (pi the angle between 
the impact parameter b± and dipole r*j. To obtain the above result we have used 
the longitudinal momentum fractions Z\ = z 2 = 0.5. These values, however, do not 
affect the asymptotic behaviour of R(s). 

With the above energy and impact parameter dependence of the x P -function, 
the black disc radius R(s) can now be obtained as follows: For cm. energies 
y/s 3> 10 6 GeV, perturbative correlations \ P dominate the dipole-dipole scattering 
amplitude S^d hi (4.2.1) because of the stronger energy dependence (e p = 0.73 > 
e = 0.125) as compared with non-perturbative correlations x ■ Consequently, 
cos(cx P ) oscillates much more than cos(cx ) in the integration over the dipole sizes 
and orientations (c = 1/3 or 2/3), i.e., only cos(cx P ) needs to be considered here. 
The condition for the black disc radius R(s) results from the fact that for \b\ < R(s) 
the functions cos(cx P ) have to average to zero when integrated over the dipole ori- 
entations in the profile function (4.2.1). This is the case when the main contribution 
to the x P -function (4.2.4) at large impact parameters fulfilles the relation 



Taking the logarithm on both sides of Eq (4.2.5), one obtaines for the leading con- 
tribution to the black disc radius R(s) at ultra-high energies 



where So denotes the reference energy, mg the gluon mass, and Mq B = Ivtlq the 
"glueball mass" obtained from our perturbative contribution [17]. The logarithmic 
increase of R(s) is equivalent to the Froissart bound [7] since oj$(s) = 2irR(s) 2 at 
asymptotic energies (see the following chapter). 

p 

For the derivation of the Froissart bound the powerlike energy dependence s e as 
well as the exponential fall off exp(— mg b) are crucial. This provides an a posteriori 
justification for the powerlike energy ansatz used in this work. The exponential 
- decrease exists because of the non-zero gluon mass tuq which has been introduced 
in the gluon propagator as a cutoff to suppress the effect of perturbative physics in 
the non-perturbative region of small momenta. 

The coefficient e p / (2mc) of the ln(s) term contains the energy exponent e p of the 
perturbative contribution and the semiperturbative parameter mc which determines 
the interplay between perturbative and non-perturbative physics. This coefficient 




(4.2.5) 




(4.2.6) 
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is universal for all hadron interactions computed within our model. Different coef- 
ficients have been used in the literature: Lukaszuk and Martin [102] have used the 
factor \pKjm v (m n is the pion mass), Heisenberg [103] has obtained by considering 
the dynamics of the scattering A /7r/(2m 7r ), Ferreiro et al. [104] have deduced from 
the Color Glass Condensate model the result a/tt/2 (4 \n(2)a s N c ) / (irm^) with the 
BFKL-exponent 4 ln(2)a s iV c /7r, and Dosch et al. [105] have updated Heisenberg's 
approach to get y/n / \2Mgb)- So far the coefficient of the ln(s) term is not conclu- 
sive. 

A universal black disc radius R(s) means that all hadrons are of the same size at 
asymptotic energies. This behavior starts becoming visible in the profile functions 
shown in Fig. 4.1 at the highest energies. The universal increase exhibits itself most 
explicitly in total hadronic cross sections which we discuss extensively in chapter 5. 



4.3 Profile Function for Photon-Proton Scattering 

In this section we compute the photon-proton scattering amplitude which is very 
appropriate to study the dependence of profile functions on the size of the scattered 
particles by varying the photon virtuality. 

The profile function for a longitudinal photon scattering off a proton p 
J y * L p(s,\b±\,Q 2 ) = [dz 1 d 2 r 1 Jdz 2 d 2 r2\^i(z 1 ,r 1 ,Q 2 )\ 2 \^ p (z2,r 2 )\ 2 

- S DD {b ± , s,z 1 ,r 1 ,z 2 ,r 2 ) (4.3.1) 

is calculated with the longitudinal photon wave function \xjj^*(zi, fi, Q 2 )\ 2 given 
in (B.2.2). In this way, the profile function (4.3.1) is ideally suited for the investiga- 
tion of dipole size effects since the photon virtuality Q 2 determines the transverse 
size of the dipole (3.3.1) into which the photon fluctuates before it interacts with 
the proton. 

Figure 4.2 shows the \b±\ dependence of the profile function J 7 £ p (s, \b±\,Q 2 ) di- 
vided by a/n for cm. energies from 10 GeV to 10 9 GeV and photon virtualities 
Q 2 = 1, 10, and 100 GeV 2 , where a is the fine-structure constant. One clearly sees 
that the qualitative behavior of this rescaled profile function is similar to the one for 
hadron-hadron scattering. However, the black disc limit induced by the underlying 
dipole-dipole scattering depends on the photon virtuality Q 2 and is given by the 
normalization of the longitudinal photon wave function 

J™p(Q 2 )= fdzd 2 r\^ L {z^Q 2 )\ 2 (4.3.2) 
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Figure 4.2: The profile function for a longitudinal photon scattering off a proton 
J 7 £p(s, Q 2 ) divided by a/n is shown versus the impact parameter for photon 
virtualities Q 2 = 1, 10, and 100 GeV 2 and cm. energies y/s from 10 GeV to 10 9 GeV. The 
value of the black disc limit J^l^iQ 2 ) and the width of the profiles depend on Q 2 . 
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since the proton wave function is normalized to one. 



The photon virtuality Q 2 does not only determine the absolute value of the 
black disc limit and the cm. energy at which it is reached but also the width of 
the profiles at fixed \fs as illustrated in Fig. 4.2. With increasing resolution Q 2 , 



i.e. decreasing dipole sizes, 



711 



oc l/Q , the absolute value of the black disc 



limit grows and higher energies are needed to reach this limit. This can be seen 
most explicitly in Fig. 4.3 where the y/s dependence of J 7 £ P (s, = 0,Q 2 ) divided 
by ajix is presented for Q 2 = 1, 10, andlOOGeV 2 . The growth of the absolute 
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Figure 4.3: The profile function for a longitudinal photon scattering off a proton 
J 7 * p (s, Q 2 ) divided by a/ir is shown versus the cm. energy yfs at zero impact param- 
eter = 0) for photon virtualities Q 2 = 1, 10, and 100 GeV 2 . 



value of the black disc limit is simply due to the normalization of the longitudinal 
photon wave function while the requirement of higher energies to reach this limit is 
due to the decreasing interaction strength with decreasing dipole size. The latter 
explains also why the energies needed to reach the black disc limit in tt p and K p 

2 Note that the Bjorken x at which the black disc limit is reached decreases with increasing 
photon virtuality Q 2 . (See also Fig. 4.4) 
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scattering are slightly higher than in pp (pp) scattering. In comparison to the hadron- 
hadron profile function, Fig. 4.2 shows explicitly that the width of the 7£p profiles 
at energies somewhat beyond the black disc limit decreases with decreasing photon 
size |r 7 *| 2 oc l/Q 2 - Furthermore, comparing j^P scattering at Q 2 = 1 GeV 2 with pp 
scattering quantitatively, the black disc limit JI^ p x (Q 2 = 1 GeV 2 ) = 0.00164 is about 
three orders of magnitude smaller because of the photon wave function normalization 
(oc a/71"). At |6j_| =0 it is reached at an energy of ^/s ~ 10 8 GeV, which is about 
two orders of magnitude higher because of the smaller dipoles involved. 

The way in which the profile function J 7 * p (s, Q 2 ) approaches the black disc 
limit at high energies depends on the shape of the proton and longitudinal photon 
wave function at small dipole sizes | ^1,2 1 - At high energies, dipoles of typical sizes 
< |ri >2 | < Ro (so/ s) 1 / 4 give the main contribution to S^* p = 1 — J 7 £ p because 
of (2.6.4) and the fact that the contribution of the large dipole sizes averages to zero 
upon integration over the dipole orientations, cf. also (4.2.2). Since S^* p is a measure 
of the proton transmittance, this means that only small dipoles can penetrate the 
proton at high energies. Increasing the energy further, even these small dipoles are 
absorbed and the black disc limit is reached. However, the dependence of the profile 
function on the short distance behavior of normalizable wave functions is weak which 
can be understood as follows. Because of color transparency, the eikonal functions 
X NP (s) and x P ( s ) are small for small dipole sizes < jri^l < -Ro (so/s) 1 ^ at large 
yfs. Consequently, Sdd ~ 1 and 

Jr L p( s ' \b±\,Q 2 ) 

1 r c (s) 1 r c (s) 

JT p x (Q 2 )-^ 2 [d Zl fdnnlMzi^Q 2 )] 2 [dz 2 [dr 2 r 2 \^ p (z 2j r 2 )\ 2 (4.3.3) 





where r c (s) ~ i?o (so/s) 1 ^. Clearly, the linear behavior from the phase space fac- 
tors ri t 2 dominates over the r 12 -dependence of normalizable wave functions. 3 More 
generally, for any profile function involving normalizable wave functions, the way in 
which the black disc limit is approached depends only weakly on the short distance 
behavior of the wave functions. 



3 For our choice of the wave functions in (4.3.3), one sees very explicitly that the specific Gaus- 
sian behavior of \ifj p (z2, ^2)| 2 and the logarithmic short distance behavior of \ip-y* (zi, 7*1, Q )| is 
dominated by the phase space factors r\^- 
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4.4 A Scenario for Gluon Saturation 



In this section we estimate the impact parameter dependent gluon distribution of 
the proton xG(x, Q 2 , \b±\). Using a leading twist, next-to-leading order QCD rela- 
tion between xG(x, Q 2 ) and the longitudinal structure function F L (x, Q 2 ), we relate 
xG(x,Q 2 , \b±\) to the profile function J 7 £ P (s = Q 2 /x, \b±\,Q 2 ) and find low- a; sat- 
uration of xG(x,Q 2 , \b±\) as a manifestation of S-matrix unitarity. The resulting 
xG(x, Q 2 , \b±\) is, of course, only an estimate since our profile function contains also 
higher twist contributions. Furthermore, in the considered low-x region, the lead- 
ing twist, next-to-leading order QCD formula may be inadequate as higher twist 
contributions [106] and higher order QCD corrections [107, 108] are expected to be- 
come important. Nevertheless, still assuming a close relation between Fl(x, Q 2 ) and 
xG(x, Q 2 ) at low x, we think that our approach provides some insight into the gluon 
distribution as a function of the impact parameter and into its saturation. 

The gluon distribution of the proton xG(x, Q 2 ) has the following meaning: 
xG(x, Q 2 )dx gives the momentum fraction of the proton which is carried by the glu- 
ons in the interval [x, x+dx] as seen by probes of virtuality Q 2 . The impact parameter 
dependent gluon distribution xG(x, Q 2 , \b±\) is the gluon distribution xG(x, Q 2 ) at a 
given impact parameter \b±\ so that 

xG(x,Q 2 ) = fd 2 b ± xG(x,Q 2 ,\b ± \) . (4.4.1) 



In leading twist, next-to-leading order QCD, the gluon distribution xG(x, Q 2 ) is 
related to the structure functions F L (x, Q 2 ) and F 2 (x, Q 2 ) of the proton [109] 



F L (x,Q' 



7T 



1 dy I x 

y \y 



r J X 



1 dy I x 

y \y 



yG(y,Q* 



(4.4.2) 

where Y2f e / is a flavor sum over the quark charges squared. For four active flavors 
and x < 10~ 3 , (4.4.2) can be approximated as follows [110] 



xG(x, 



5.8 



3vr 
4a, 



F L (0.417x,Q s 



1.97 



F 2 (0.75x,Q 2 



(4.4.3) 



For typical A QCD = 100 - 300 MeV and Q 2 = 50 - 100 GeV 2 , the coefficient of 
Fl in (4.4.3), 37r/(4a s ) = O(10), is large compared to the one of F 2 . Taking into 
account also the values of F 2 and Fl, in this Q 2 region and for x < 10~ 3 , the gluon 
distribution is mainly determined by the longitudinal structure function. The latter 
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can be expressed in terms of the profile function for longitudinal photon-proton 
scattering using the optical theorem (cf. (5.1.1)) 

F L (x,Q 2 ) = ^a t °Ux,Q 2 ) = ^ 2 [d 2 b ± J rLP (x, \b ± \,Q 2 ) , (4-4.4) 

where the s-dependence of the profile function is rewritten in terms of the Bjorken 
scaling variable, x = Q 2 /s. Neglecting the F2 term in (4.4.3), consequently, the 
gluon distribution reduces to 

xG(x, Q 2 ) « 1.305 - [d 2 b ± J 7 * p (0.417x, \b ± \,Q 2 ) . (4.4.5) 
n 2 a s a J L 

Comparing (4.4.1) with (4.4.5), it seems natural to relate the integrand of (4.4.5) to 
the impact parameter dependent gluon distribution 

xG(x, Q 2 , |6 ± |) « 1.305 ^- - J 7 * p (0.417x, \b ± \, Q 2 ) . (4.4.6) 



The black disc limit of the profile function for longitudinal photon-proton scat- 
tering (4.3.2) imposes accordingly an upper bound on xG(x,Q 2 , \b±\) 

xG(x,Q 2 ,\b ± \) < xG max (Q 2 ) w 1.305-^- J"7(Q 2 ) , (4.4.7) 

which is the low-x saturation value of the gluon distribution xG(x,Q 2 , \b±\) in our 
approach. With it J™®® (Q 2 ) / a « 1 as shown in Fig. 4.3, a compact approximation 
of (4.4.7) is obtained 

xG(x,Q 2 , \b±\) < xG max {Q 2 ) w , (4.4.8) 

which is consistent with the results in [108, 111, 112] and indicates strong color-field 
strengths G a ~ 1/ ^/oT s as well. 

According to our relations (4.4.6) and (4.4.7), the blackness described by the pro- 
file function is a measure for the gluon density and the black disc limit corresponds to 
the maximum gluon density reached at the impact parameter under consideration. 
In accordance with the behavior of the profile function J 7 £ P , see Fig. 4.2, the gluon 
distribution xG(x, Q 2 , \b±_\) decreases with increasing impact parameter for given 
values of x and Q 2 . The gluon density, consequently, has its maximum in the geo- 
metrical center of the proton, i.e. at zero impact parameter, and decreases towards 
the periphery. With decreasing x at given Q 2 , the gluon distribution xG(x, Q 2 , \b±\) 
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increases and extends towards larger impact parameters just as the profile function 
J 7 £ P for increasing s. The saturation of the gluon distribution xG(x,Q 2 , \b±\) sets 

in first in the center of the proton (\b±\ =0) at very small Bjorken x. 

In Fig. 4.4, the gluon distribution xG(x, Q 2 , \b±\ = 0) is shown as a function 
of x for Q 2 = 1, 10, andlOOGeV 2 , where the relation (4.4.6) has been used also 
for low photon virtualities. Evidently, the gluon distribution xG(x, Q 2 , \b±\ = 0) 
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Figure 4.4: The gluon distribution of the proton at zero impact parameter 
xG(x, Q 2 , \b±\ = 0) is shown as a function of x for Q 2 = 1, 10, and 100 GeV 2 . The re- 
sults are obtained within the approximation (4.4.6). 

saturates at very low values of x < 10~ 10 for Q 2 > 1 GeV 2 . The photon virtuality 
Q 2 determines the saturation value (4.4.7) and the Bjorken-x at which it is reached 
(cf. also Fig. 4.2). For larger Q 2 , the low-x saturation value is larger and is reached at 
smaller values of x, as claimed also in [113]. Moreover, the growth of xG(x, Q 2 , \b±\ = 
0) with decreasing x becomes stronger with increasing Q 2 . This results from the 
stronger energy increase of the perturbative component, e p = 0.73, that becomes 
more important with decreasing dipole size. 

According to our approach, the onset of the xG(x,Q 2 , |6jJ)-saturation appears 
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for Q 2 > 1 GeV 2 at x < 10~ 10 , which is far below the x-region accessible at HERA 
(x > 10~ 6 ). Even for THERA (x > 10~ 7 ), gluon saturation is not predicted for Q 2 > 
1 GeV 2 . Munier et al. have also shown in a recent publication [114] that the black 
disc limit is not reached in the energy range of HERA. However, since the HERA 
data can be described by models with and without saturation embedded [113], the 
present situation is not conclusive. 

Note that the 5*- matrix unitarity condition (4.1.3) together with (4.4.6) requires 
the saturation of the impact parameter dependent gluon distribution xG(x, Q 2 , \b±\) 
but not the saturation of the integrated gluon distribution xG(x, Q 2 ). Due to mul- 
tiple gluonic interactions in our model, this requirement is fulfilled, as can be seen 
from Fig. 4.2 and relation (4.4.6). Indeed, approximating the gluon distribution 
xG(x, Q 2 , \b±\) in the saturation regime of very low x by a step-function 

xG(x,Q\\b i _\)^xG max (Q 2 )Q(R(x,Q 2 )-\b 1 _\) , (4.4.9) 

where R(x, Q 2 ) denotes the black disc radius (or the full width at half maximum 
of the profile function), one obtains with (4.4.1), (4.4.7) and (4.4.8) the integrated 
gluon distribution 

xG(x,Q 2 ) « 1.305 g2jR2(x ' g2) -J— (Q 2 ) « Q2jR2(X,g2) , (4.4.10) 

which does not saturate because of the increase of the effective proton radius R(x, Q 2 ) 
with decreasing x. Nevertheless, although xG(x, Q 2 ) does not saturate, the satu- 
ration of xG(x, Q 2 , \bj_\) leads to a slow-down from powerlike to squard logarithmic 
(see (4.2.6)) growth of xG(x, Q 2 ) towards small x. 5 Interestingly, our result (4.4.10) 
coincides with the result of Mueller and Qiu [108]. 

Finally, it must be emphasized that the low-x saturation of xG(x,Q 2 , \b±\), re- 
quired in our approach by the S'-matrix unitarity, is realized by multiple gluonic in- 
teractions. In other approaches that describe the evolution of the gluon distribution 
with varying x and Q 2 , gluon recombination leads to gluon saturation [107, 108, 116— 
118], which is reached when the probability of a gluon splitting up into two is equal 
to the probability of two gluons fusing into one. A more phenomenological under- 
standing of saturation is attempted in [71, 119]. 



4 So far, the most striking hint for saturation in the present HERA data at x ~ 10~ 4 and Q 2 < 
2 GeV 2 has been the turnover of dF^{x, Q 2 ) / dln(Q 2 ) towards small x in the Caldwell plot [115], 
which is still a controversial issue due to the correlation of Q 2 and x values. 

5 This is analogous to the slow down of the total pp cross section at high cm. energy as soon as 
the corresponding profile function J pp (s, \b±\) reaches its black disc limit as shown in Sec. 5.1. 
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Chapter 5 



Comparison with Data and 
Saturation Effects in Observables 



In this chapter we show the phenomenological performance of our model. We com- 
pute total, differential, and elastic cross sections, structure functions, and diffractive 
slopes for hadron-hadron, photon-proton, and photon-photon scattering, compare 
the results with experimental data including cosmic ray data, and provide predic- 
tions for future experiments. Making use of the explicit saturation of the impact 
parameter profiles at the black disc limit, the corresponding energy values, and the 
logarithmic rise of the black disc radius studied in the preceding chapter, we show 
explicitly manifestations of S'-matrix unitarity limits in experimental observables. 

Using the T-matrix element given in (2.6.6), we compute the pomeron contribu- 
tion to pp, pp, 7r ± jo, K ± p, 7*p, and 77 reactions in terms of the universal dipole- 
dipole scattering amplitude Sdd- This allows one to compare reactions induced by 
hadrons and photons in a systematic way. In fact, it is our aim to provide a unified 
description of all these reactions and to show in this way that the pomeron contribu- 
tion to the above reactions is universal and can be traced back to the dipole-dipole 
scattering amplitude Sdd- 

Our model describes pomeron (C = +1 gluon exchange) but neither odderon 
(C = — 1 gluon exchange) nor reggeon exchange (quark-antiquark exchange) as 
discussed in Sec. 2.5. Only in the computation of the hadronic total cross sections 
the reggeon contribution is added [29,120]. This improves the agreement with the 
data for y/s < 100 GeV and describes exactly the differences between hadron-hadron 
and antihadron-hadron reactions. 

The model parameters have been adjusted in fits to the high-energy scattering 
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data shown in this chapter. The resulting parameter set given in Sec. 2.7 and 
Appendix B is used throughout this work. 



5.1 Total Cross Sections 



The total cross section for the high-energy reaction ah — > X is related via the optical 
theorem to the imaginary part of the forward elastic scattering amplitude and can 
also be expressed in terms of the profile function (4.1.2) 

o%{a) = - s lmT(s,t = 0) = 2 J ' d 2 b ± J ab (s,\b ± \) , (5.1.1) 

where a and b label the initial particles whose masses were neglected as they are 
small in comparison to the cm. energy y'i. 

We compute the pomeron contribution to the total cross section, c*^ ,JP (s), from 
the T-matrix (2.6.6), as explained above, and add only here a reggeon contribution 
of the form [29, 120] 

i + id / S \ -0.4525 

tot.IRr \ -\r I 



where X ab depends on the reaction considered: X pp = 56.08 mb, X pp = 98.39 mb, 
X n + P = 27.56 mb, X n - p = 36.02 mb, X K + P = 8.15 mb, X K - p = 26.36 mb, X lp = 
0.129 mb, and X 77 = 605 nb. Accordingly, we obtain the total cross section 

<#(«)= <#^(«) + ^(*) (5-1-3) 
for pp, pp, 7r ± jo, K ± p, 7p and 77 scattering. 

The good agreement of the computed total cross sections with the experimental 
data is shown in Fig. 5.1. Here the solid lines represent the theoretical results for 
pp, n + p, K + p, 7p, and 77 scattering and the dashed lines the ones for pp, ix~p, 
and K~p scattering. The dot-dashed line represents 27r(e p /2mc) 2 ln 2 (s/s ) with 
y/1 = 20 GeV. The pp, pp, n^p, K ± p, jp [3] and 77 data [121] taken at accelerators 
are indicated by the closed circles while the closed squares (Fly's eye data) [122] 
and the open circles (Akeno data) [123] indicate cosmic ray data. Concerning the 
photon-induced reactions, only real photons are considered which are, of course, 
only transverse polarized. 



The prediction for the total pp cross section at LHC (^/s = 14TeV) is a, 



tot 

pp 



114.2 mb in good agreement with the cosmic ray data. Compared with other works, 
our LHC prediction is close to the one of Block et al. [124], a p p = 108 ± 3.4 mb, but 
considerably larger than the one of Donnachie and Landshoff [29], o 1 ^ = 101.5 mb. 
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Figure 5.1: The total cross section a tot is shown as a function of the cm. energy yfs for 
pp, pp, tt^p, K^p, jp and 77 scattering. The solid lines represent the model results for 
pp, 7r + p, K + p, 7p and 77 scattering and the dashed lines the ones for pp, ir~p, and K~p 
scattering. The dot-dashed line represents 27r(e p /2ma) 2 ln 2 (s/so) with \/I = 20 GeV. 
The pp, pp, n^p, K p, 7p [3] and 77 data [121] taken at accelerators are indicated by the 
closed circles while the closed squares (Fly's eye data) [122] and the open circles (Akeno 
data) [123] indicate cosmic ray data. The arrows at the top point to the LHC energy, 
sfs = 14TeV, and to the onset of the black disc limit in pp (pp) reactions, sfs ~ 10 6 GeV. 
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The differences between ab and ab reactions for y/s < 100 GeV result solely from 
the different reggeon contributions which die out rapidly as the energy increases. The 
pomeron contribution to ab and ab reactions is, in contrast, identical and increases 
as the energy increases. It thus governs the total cross sections for y/s > 100 GeV 
where the results for ab and ab reactions coincide. 

The differences between pp (pp) , n^p, and K ± p scattering result from the differ- 
ent transverse extension parameters, S p = 0.86 fm > S n = 0.607 fm > Sk = 0.55 fm, 
cf. Appendix B. Since a smaller transverse extension parameter favors smaller 
dipoles, the total cross section becomes smaller, and the short distance physics 
described by the perturbative component becomes more important and leads to a 
stronger energy growth due to e p = 0.73 > e NP = 0.125. In fact, the ratios <r^f/<7*p 
and ct^/ctkj, converge slowly towards unity with increasing energy as can already 
be seen in Fig. 5.1. 

For real photons, the transverse size is governed by the constituent quark masses 
iTif(Q 2 = 0), cf. Appendix B, where the light quarks have the strongest effect, i.e. 
cr*p oc l/m^ d and <r*°' oc l/m^ d . Furthermore, in comparison with hadron-hadron 
scattering, there is the additional suppression factor of a for 7p and a 2 for 77 
scattering coming from the photon-dipole transition. In the 77 reaction, also the 
box diagram contributes [120, 125] but is neglected since its contribution to the total 
cross section is less than 1% [74]. 

It is worthwhile mentioning that total cross sections for pp (pp), 7r p, and K ± p 
scattering do not depend on the width Azh of the longitudinal quark momentum 
distribution in the hadrons (B.l.l) since the underlying dipole-dipole scattering 
process is independent of the longitudinal quark momentum fraction for t = 0. 
On the two-gluon-exchange level, this is shown analytically in chapter 3. 

Saturation effects as a manifestation of the S"-matrix unitarity can be seen in 
Fig. 5.1. Having investigated the profile function for hadron-hadron scattering, we 
know that this profile function becomes higher and broader with increasing energy 
until it saturates the black disc limit first for zero impact parameter = 0) at 
\/s pd 10 6 GeV. Beyond this energy, the profile function cannot become higher but 
expands towards larger values of Consequently, the total cross section (5.1.1) 
increases no longer due to the growing blackness at the center but only due to the 
transverse expansion of the hadrons. This tames the growth of the total hadronic 
cross sections as can be seen for cm. energies beyond \fs ~ 10 6 GeV in Fig. 5.1. 

At energies far beyond the onset of the black disc limit at zero impact parameter, 
the profile function can be approximated by 

JT OX ^ |5lI) = N a N b Q (R(s) - \b ± \) (5.1.4) 
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where N a ^ denotes the normalization of the wave functions of the scattered particles 
and R(s) the black disc radius defined in Sec. 4.2 that reflects the effective radii of the 
interacting particles. Thus, the energy dependence of the total cross section (5.1.1) 
is driven exclusively by the increase of the transverse extension of the particles R(s) 



which is known as geometrical scaling [126,127]. Introducing the analytical result 
for the black disc radius R(s) at asymptotic energies (4.2.6) in the above equation, 
one obtaines the ln 2 -growth 



which coincides with the Froissart bound [7]. For hadron wave function normaliza- 
tions N a)b = 1, the total cross section (5.1.6) is independent of the hadron species 
involved in the hadron-hadron scattering at asymptotic energies. The hadronic cross 
sections start joining already at the highest energies shown in Fig. 5.1. Also for pho- 
tons of different virtuality Q\ and Q\ one can check that the ratio of the total cross 
sections (Qf) / 'cr^ p (Ql) converges to unity at asymptotic energies in agreement 
with the conclusion in [128]. 

The total cross section (5.1.6) obtained with the values for e p and vtlq from 
Sec. 2.7, N a> b = 1, and the cm. energy y/Io = 20 GeV is shown as a dot-dashed 
line for ^/s > 600 GeV in Fig. 5.1. Evidently, the so obtained logarithmic growth 
coincides with the numerical result for the pp cross section at high energies. In fact, 
a transition from a powerlike to an ln 2 -increase of total cross sections seems to set in 
at about \/s = 10 6 GeV as visible in Fig. 5.1. Following the recent publications [105, 
129], already fits to available forward scattering data support the universal In (s)- 
dependence of the cross sections. 



5.2 Proton Structure Function 

The total cross section for the scattering of a transverse (T) and longitudinally 
(L) polarized photon off the proton, o l °* p (x, Q 2 ), at photon virtuality Q 2 and cm. 

energy 1 squared, s = Q 2 /x, is equivalent to the structure functions of the proton 



a%{s) = 2nN a N b R(s) 2 , 



(5.1.5) 




(5.1.6) 




(x,Q 2 ) 



(5.2.1) 



Here, ^fs refers to the cm. energy in the j*p system. 
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and 



F 2 (x,Q 2 ) = F T (x,Q 2 ) + F L (x,Q 2 ) . 



(5.2.2) 



Reactions induced by virtual photons are particularly interesting because the 
transverse separation of the quark-antiquark pair that emerges from the virtual 
photon decreases as the photon virtuality increases, cf. Eq. (3.3.1). With increasing 
virtuality, one probes therefore smaller transverse distance scales of the proton. 

In Fig. 5.2, the Q 2 -dependence of the total 7*p cross section 



is presented, where the model results (solid lines) are compared with the experi- 
mental data for cm. energies from yfs = 20 GeV up to y/s = 245 GeV. Note the 
indicated scaling factors at different y/s values. The low energy data at y/s = 20 GeV 
are from [130] while the data at higher energies have been measured at HERA by 
the HI [131] and ZEUS collaboration [132]. At Q 2 = 0.012 GeV 2 , also the photo- 
production (Q 2 = 0) data from [133] are displayed. 

In the window shown in Fig. 5.2, the model results are in reasonable agreement 
with the experimental data. The total 7*p cross section levels off towards small val- 
ues of Q 2 as soon as the photon size |r 7 |, i.e the resolution scale, becomes comparable 
to the proton size. Our model reproduces this behavior by using the perturbative 
photon wave functions with Q 2 -dependent quark masses, mf(Q 2 ), that interpolate 
between the current (large Q 2 ) and the constituent (small Q 2 ) quark masses as ex- 
plained in detail in Appendix B. The decrease of a l °* p with increasing Q 2 results 
from the decreasing dipole sizes since small dipoles do not interact as strongly as 
large dipoles. 

The x-dependence of the computed proton structure function F 2 (x, Q 2 ) is shown 
in Fig. 5.3 for Q 2 = 0.3, 2.5, 12 and 120 GeV 2 in comparison to the data measured by 
the HI [134] and ZEUS [135] detector. Within our model, the increase of F 2 (x,Q 2 ) 
towards small Bjorken x becomes stronger with increasing Q 2 in agreement with the 
trend in the HERA data. This behavior results from the fast energy growth of the 
perturbative component that becomes more important with increasing Q 2 due to 
the smaller dipole sizes involved. 

As can be seen in Fig. 5.3, the data show a stronger increase with decreasing x 
than the model outside the low-Q 2 region. This results from the weak energy boost 
of the non-perturbative component that dominates F 2 (x, Q 2 ) in our model. In fact, 
even for large Q 2 the non-perturbative contribution overwhelms the perturbative 
one, which explains also the overestimation of the data for x > 10~ 3 . 



a%(s,Q 2 ) = a% p (s,Q 2 ) + a% p (s,Q 2 ) 



(5.2.3) 
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Q 2 [GeV 2 ] 

Figure 5.2: The total 7*p cross section, ai°^ p {s, Q 2 ), is shown as a function of the photon 
virtuality Q 2 for cm. energies from \/s = 20 GeV to y/s = 245 GeV, where the model 
results (solid lines) and the experimental data at different y/s values are scaled with the 
indicated factors. The low energy data at yfs = 20 GeV are from [130], the data at higher 
energies from the HI [131] and ZEUS collaboration [132]. The photoproduction (Q 2 = 0) 
data from [133] are displayed at Q 2 = 0.012 GeV 2 . 
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Figure 5.3: The x-dependence of the computed proton structure function F 2 {x, Q 2 ) (solid 
line) is shown for Q 2 = 0.3, 2.5, 12 and 120 GeV 2 in comparison to the data measured by 
the HI [134] and ZEUS [135] detector, and the low energy data at </s = 20 GeV from [130]. 



This problem is typical for the SVM model applied to the scattering of a small 
size dipole off a proton. In an earlier application by Riiter [21], an additional cut-off 
was introduced to switch from the non-perturbative to the perturbative contribution 
as soon as one of the dipoles is smaller than r cut = 0.16 fm. This yields a better 
agreement with the data also for large Q 2 but leads to a discontinuous dipole- 
proton cross section. In the model of Donnachie and Dosch [74], a similar SVM- 
based component is used also for dipoles smaller than R c = 0.22 fm with a strong 
energy boost instead of a perturbative component. Furthermore, their SVM-based 
component is tamed for large Q 2 by an additional a s (Q 2 ) factor. 

We did not follow these lines in order to keep a continuous, Q 2 -independent 
dipole-proton cross section and, therefore, cannot improve the agreement with the 
F 2 (x, Q 2 ) data without losing quality in the description of hadronic observables. 
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Since our non-perturbative component relies on lattice QCD, we are more confi- 
dent in describing non-perturbative physics and, thus, put more emphasis on the 
hadronic observables. Admittedly, our model misses details of the proton structure 
that become visible with increasing Q 2 . In comparison, most other existing models 
provide neither the profile functions nor a simultaneous description of hadronic and 
7*-induced processes. For example, the phenomenological model of Golec-Biernat 
and Wiisthoff [71, 72] provides a successful and economical description of the 7*p 
reactions but cannot be applied to hadron-hadron reactions. 



5.3 Slope Parameter B of Elastic Forward Scattering 



The local slope of elastic scattering B(s,t) is defined as 

~da el 



dt 



(s,t) 



(5.3.1) 



and, thus, characterizes the diffractive peak of the differential elastic cross section 
do~ el /dt(s,t) discussed below. Here, we concentrate on the slope for elastic forward 
(t = 0) scattering also called slope parameter 

B(.) := B(s, t = 0) = llf^^lM = V) , (5.3.2) 

2 jcn ± J(s,\b±\) 2 

which measures the rms interaction radius (b 2 ) of the scattered particles, and does 
not depend on the opacity. 

We compute the slope parameter with the profile function from the T-matrix 
(2.6.6) and neglect the reggeon contributions, which are relevant only at small cm. 
energies, so that the same result is obtained for ab and db scattering. 

In Fig. 5.4, the resulting slope parameter B(s) is shown as a function of y/s for 
pp and pp scattering (solid line) and compared with the pp (open circles) and pp 
(closed circles) data from [136-138]. As expected from the opacity independence of 
the slope parameter (5.3.2), saturation effects as seen in the total cross sections do 
not occur. Indeed, one observes an approximate B(s) oc R 2 (s) cx ln 2 (s) growth for 
x/s > 10 4 GeV. This behavior agrees, of course, with the transverse expansion of 
Jpp(s, \b±\) for increasing yfs shown in Fig. 4.1. Analogous results are obtained also 
for np and Kp scattering. 

For the good agreement of our model with the data, the finite width of the lon- 
gitudinal quark momentum distribution in the hadrons, i.e. Az p , Az n , and Azk ^ 
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in (B.l.l), is important as the numerator in (5.3.2) depends on this width. In 
fact, B(s) comes out more than 10% smaller with Az p , Az w , andA^ = 0. The 
A^-dependence of the non-forward scattering amplitude computed analytically on 
the two-gluon exchange level is shown explicitly in Appendix D. Furthermore, a 
strong growth of the perturbative component, e p = 0.73, is important to achieve 
the increase of B(s) for y^s > 500 GeV indicated by the data. 
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Figure 5.4: The elastic slope parameter B(s) is shown as a function of the cm. 
energy yfs for pp and pp forward (t = 0) scattering. The solid line represents the 
model result that is compared with the data for pp (open circles) and pp (closed 
circles) reactions from [136-138]. 

It must be emphasized that only the simultaneous fit of the total cross section and 
the slope parameter provides the correct shape of the profile function J(s, \b±\). This 
shape leads then automatically to a good description of the differential elastic cross 
section da el /dt(s,t) as demonstrated below. Astonishingly, only few phenomeno- 
logical models provide a satisfactory description of both quantities [70, 124]. In the 
approach of [14] , for example, the total cross section is described correctly while the 
slope parameter exceeds the data by more than 20% already at y/s = 23.5 GeV and, 
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thus, indicates deficiencies in the form of J(s, \b±\). 



5.4 Differential Elastic Cross Sections 



The differential elastic cross section obtained from the squared absolute value of the 
T-matrix element 

^{s,t) = —AT^M (5.4.1) 
at lbirs z 

can be expressed for our purely imaginary T-matrix (2.6.6) in terms of the profile 
function 



da el , , 1 



dt ^ t) = T, 



Jd 2 b±e icL ^ J(s,\b ± \) 



(5.4.2) 



and is, thus, very sensitive to the transverse extension and opacity of the scattered 
particles. Equation (5.4.2) reminds of optical diffraction, where J(s, \b±\) describes 
the distribution of an absorber that causes the diffraction pattern observed for inci- 
dent plane waves. 

In Fig. 5.5, the differential elastic cross section computed for pp and pp scattering 
(solid line) is shown as a function of \t\ = at y/s = 23.5, 30.7, 44.7, 63, 546, and 
1800 GeV and compared with experimental data (open circles), where the pp data 
at ■y/s = 23.5, 30.7, 44.7, and 63 GeV were measured at the CERN ISR [126], the pp 
data at y/s = 546 GeV at the CERN SppS [137], and the pp data at y/s = 1.8 TeV at 
the Fermilab Tevatron [138, 139]. The prediction of our model for the pp differential 
elastic cross section at the CERN LHC, y/s = 14 TeV, is given in Fig. 5.6. 

For all energies, the model reproduces the experimentally observed diffraction 
pattern, i.e the characteristic diffraction peak at small \t\ and the dip structure 
at medium \t\. As the energy increases, also the shrinking of the diffraction peak is 
described which reflects the rise of the slope parameter B(s,t = 0) already discussed 
above. The shrinking of the diffraction peak comes along with a dip structure that 
moves towards smaller values of \t\ as the energy increases. This motion of the dip 
is also described approximately. 

The dip in the theoretical curves reflects a change of sign in the T-matrix el- 
ement (2.6.6). As the latter is purely imaginary, it is not surprising that there 
are deviations from the data in the dip region. Here, the real part is expected to 
be important [139] which is in the small |t| region negligible in comparison to the 
imaginary part. 

The difference between the pp and pp data, a deep dip for pp but only a bump 
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Figure 5.5: The differential elastic cross section for pp and pp scattering is shown as 
a function of \t\ up to 2.5 GeV 2 . The result of our model, indicated by the solid line, 
is compared for ^ = 23.5, 30.7, 44.7, and 63 GeV with the CERN ISR pp data [126], 
for = 546 GeV with the CERN SppS pp data [137], and for = 1.8 TeV with the 
Fermilab Tevatron pp data [138, 139], all indicated by the open circles with error bars. 
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Figure 5.6: The prediction of our model for the pp differential elastic cross section at 
LHC (y/s = 14 TeV) is shown as a function of momentum transfer \t\ up to 1 GeV 2 . 



or shoulder for pp reactions, requires a C = — 1 contribution. Besides the reggeon 
contribution at small energies, 2 one expects here an additional perturbative C = —1 
contribution such as three-gluon exchange [140,141] or an odderon [142-144]. In 
fact, allowing a finite size diquark in the (anti-)proton an odderon appears that 
supports the dip in pp but leads to the shoulder in pp reactions [144]. 

For the differential elastic cross section at the LHC energy, y/s = 14 TeV, the 
above findings suggest an accurate prediction in the small- \t\ region but a dip at a 
position smaller than the predicted value at \t\ ~ 0.35 GeV 2 . Our confidence in the 
validity of the model at small |t| is supported additionally by the total cross section 
that fixes da el /dt(s,t = 0) and is in agreement with the cosmic ray data shown 
in Fig. 5.1. Concerning our prediction for the dip position, it is close to the value 
|*| « 0.41 GeV 2 of [124] but significantly below the value \t\ w 0.55 GeV 2 of [14]. 
Beyond the dip position, the height of the computed shoulder is always above the 

2 Zooming in on the result for ^fs — 23.5 GeV, one finds further an underestimation of the data 
for all \t\ before the dip, which is correct as it leaves room for the reggeon contribution being 
non-negligible at small energies. 
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data and, thus, very likely to exceed also the LHC data. In comparison with other 
works, the height of our shoulder is similar to the one of [124] but almost one order 
of magnitude above the one of [14]. 

Considering Figs. 5.5 and 5.6 more quantitatively in the small-|t| region, one can 
use the well known parametrization of the differential elastic cross section in terms 
of the slope parameter B(s) and the curvature C(s) 

da el /dt(s, t) = da el /dt(s, t = 0) exp [B(s)t + C(s)t 2 ] . (5.4.3) 

Using B(s) from the preceding section and assuming for the moment C(s) = 0, 
one achieves a good description at small momentum transfers and energies, which 
is consistent with the approximate Gaussian shape of J pp (s, \b±\) at small energies 
shown in Fig. 4.1. The dip, of course, is generated by the deviation from the Gaus- 
sian shape at small impact parameters. According to (5.4.3), the shrinking of the 
diffraction peak with increasing energy reflects simply the increasing interaction 
radius described by B(s). 

For small energies y/s, our model reproduces the experimentally observed change 
in the slope at |t| ~ 0.25 GeV 2 [145] that is characterized by a positive curvature. 
For LHC, we find clearly a negative value for the curvature in agreement with [124] 
but in contrast to [14]. The change of sign in the curvature reflects the transition of 
J(s, \b±\) from the approximate Gaussian shape at low energies to the approximate 
step-function shape (5.1.4) at high energies. 

Important for the good agreement with the data is the longitudinal quark mo- 
mentum distribution in the proton. Besides the slope parameter, which characterizes 
the diffraction peak, also the dip position is very sensitive to the distribution width 
Az p , i.e. with Az p = the dip position appears at more than 10% lower values of |t|. 
In the earlier SVM approach [14], the reproduction of the correct dip position was 
possible without the ^-dependence of the hadronic wave functions but a deviation 
from the data in the low-|t| region had to be accepted. In this low-|i| region, we 
achieved a definite improvement with the new correlation functions (2.4.10) and the 
minimal surfaces used in our model. 

The differential elastic cross section computed for 7r ± p and K^p reactions has the 
same behavior as the one for pp (pp) reactions. The only difference comes from the 
different z-distribution widths, Az n and Azk, and the smaller extension parameters, 
S n and Sk, which shift the dip position to higher values of \t\. This is illustrated 
in Fig. 5.7, where the model results (solid line) for the n^p and K^p differential 
elastic cross section as a function of \t\ are shown at \/s = 19.5 GeV in comparison 
with experimental data (closed circles) from [146]. The deviations from the data 
towards large |i| leave room for odderon and reggeon contributions. Indeed, with a 
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Figure 5.7: The differential elastic cross section da el /dt(s,t) is shown versus the momen- 
tum transfer \t\ for 7r ± p and K^p reactions at the cm. energy = 19.5 GeV. The model 
results (solid line) are compared with the data (closed circles with error bars) from [146]. 

finite diquark size in the proton, an odderon occurs that improves the description 
of the data at large values of \t\ [147]. 



5.5 Elastic Cross Sections cr el , cr el /(7 tot , and <j tot /B 

The elastic cross section obtained by integrating the differential elastic cross section 

f-oo j el f-oo i 
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reduces for our purely imaginary T-matrix (2.6.6) to 



a el (s) = Jd 2 b ± \J(s,\b ± \)\ 2 . (5.5.2) 

Due to the squaring, it exhibits the saturation of J(s, \b±\) at the black disc limit 
more clearly than a tot (s). Even more transparent is the saturation in the following 
ratios given here for a purely imaginary T-matrix 

= If^m , (5.5.3) 



2jd% ± J(s,\b ± \) 
2fd% ± J(s,\b ± \) 



a tot 



B Jd 2 b ± \b ± \*J(s,\b ± \) 



(5.5.4) 



which are directly sensitive to the opacity of the particles. This sensitivity can be 
illustrated within the approximation 

T(s, t)=is a tot (s) exp[B(s)t/2] (5.5.5) 

that leads to the differential cross section (5.4.3) with C(s) = 0, i.e. an exponential 
decrease over \t\ with a slope B(s). As the purely imaginary T-matrix element (5.5.5) 
is equivalent to 

J(s, \b±\) = {<j tot /AnB) exp{-\b ± \ 2 /2B] = (Aa el /a tot ) eocp[— |6l| 2 /2B] , (5.5.6) 

one finds that the above ratios are a direct measure for the opacity at zero impact 
parameter 

J( s , |6 ± | = 0) = (a tot /A7rB) = (Aa el /a tot ) . (5.5.7) 

For a general purely imaginary T-matrix, T(s,t) = i s a tot g(\t\) with an arbitrary 
real-valued function g(|t|), J(s, \b±\ = 0) is given by (a el /a tot ) times a pure number 
which depends on the shape of g(|t|). 

We compute the elastic cross section a el and the ratios a el /a tot and a tot /B in 
our model without taking into account reggeons. In Fig. 5.8, the results for pp and 
pp reactions (solid lines) are compared with the experimental data (open and closed 
circles). The data for the elastic cross section are taken from [3] and the data for 
a tot and B from the references given in previous sections. For pp (pp) scattering, we 
indicate explicitly the prediction for LHC at y^i = 14 TeV and the onset of the black 
disc limit at = 10 6 GeV. The model results for irp and Kp reactions are presented 
as dashed and dotted line, respectively. For the ratios, the asymptotic limits are 
indicated: Since the maximum opacity or black disc limit governs the — > oo 
behavior, a el /a tot (a tot /B) cannot exceed 0.5 (8tt) in hadron-hadron scattering. 
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Figure 5.8: The elastic cross section a and the ratios a ja tot and a tot /B are shown as 
a function of the cm. energy ^fs. The model results for pp (pp), Tip, and Kp scattering 
are represented by the solid, dashed and dotted lines, respectively. The experimental data 
for the pp and pp reactions are indicated by the open and closed circles, respectively. The 
data for the elastic cross section are taken from [3] and the data for a tot and B from the 
references given in previous sections. 
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In the investigation of pp (pp) scattering, our theoretical curves confront suc- 
cessfully the experimental data for the elastic cross section and both ratios. At low 
energies, the data are underestimated as reggeon contributions are not taken into 
account. Again, the agreement is comparable to the one achieved in [124] and better 
than in the approach of [14], where a el comes out too small due to an underestima- 
tion of da el /dt in the low-|t| region. 

Concerning the energy dependence, a el shows a similar behavior as a tot but with a 
more pronounced flattening around y/s > 10 6 GeV. This flattening is even stronger 
for the ratios, drawn on a linear scale, and reflects very clearly the onset of the 
black disc limit. As expected from the simple approximation (5.5.7), a el /a tot and 
<j tot /B show a similar functional dependence on y/s. At the highest energy shown, 
y/s = 10 8 GeV, both ratios are still below the indicated asymptotic limits, which 
reflects that J(s, \b±\) still deviates from the step-function shape (5.1.4). The ratios 
a el /a tot and a tot / B reach their upper limits 0.5 and 8ir, respectively, at asymptotic 
energies, y/s — > oo, where the hadrons become infinitely large, completely black 
discs. 

In comparision to the pp (pp) result for a el , the irp and Kp elastic cross sections 
are smaller due to the smaller pion and kaon sizes or the corresponding narrower 
profile functions shown in Fig. 4.1. For pp (pp), irp, and Kp scattering, the ratios 
a /a tot and a tot /B representing the opacity of the interacting hadrons in their center 
(cf. Fig. 4.1 at = 0) converge for y/s > 10 6 GeV as shown in Fig. 5.8. This 
follows from the identical normalizations of the hadron wave functions that lead to 
an identical black disc limit for hadron-hadron reactions approached in the center 
of the hadrons at y/s 10 6 GeV. 
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Chapter 6 



Conclusions and Outlook 



We have developed a model that combines perturbative and non-perturbative QCD 
to compute high-energy reactions of hadrons and photons. We have investigated the 
QCD structure of non-perturbative interactions in momentum space. A decompo- 
sition of the QCD string into stringless dipoles has been found, confinement effects 
have been shown explicitly in high-energy scattering and a microscopic structure of 
unintegrated gluon distributions of hadrons and photons has been obtained. Sat- 
uration effects that manifest the unitarity of the S-matrix have been studied in 
impact parameter space of the scattering amplitude. We have diplayed explicitly 
the saturation of the impact parameter profiles for proton-proton and photon-proton 
scattering at the black disc limit at high cm. energies. The impact parameter de- 
pendent gluon distribution of the proton is found to saturate at small Bjorken x. 
We have compared the model results for pp, irp, Kp, 7*p and 77 reactions with 
experimental data and have shown saturation effects in experimental observables. 

The simplicity of minimal surfaces introduced in our model has allowed us to 
show for the first time the QCD structure of the non-perturbative dipole-dipole 
interaction in momentum space. This contribution divides into two parts: The first 
part describes the non-perturbative interaction between the quarks and antiquarks 
of the two dipoles and exhibits the structure known from perturbative two-gluon 
exchange [34,35]. The second part describes the interaction between the strings of 
the two dipoles and has a new structure originating from the geometry of the strings. 

The contribution of the confining string to the total dipole-hadron cross section 
PDh(x, \r D \) has been studied. For small dipole sizes, \f D \ — > 0, the string contri- 
bution shows color transparency, cr Dh (x, \ f D \) oc r%, as known for the perturbative 
contribution. For large dipole sizes, \f D \ > 0.5 fm, the non-perturbative contribu- 
tion increases linearly with increasing dipole size, ctd/j(x, \r D \) oc \r D \, in contrast to 
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the perturbative contribution which gives an \r D \ - independent dipole-hadron cross 
section. This linear increase is generated by the interaction of the string of the dipole 
with the hadron: the longer the string, the larger the geometric cross section with 
the hadron. String breaking, expected to stop the linear increase for \r D \ > 1 fm, 
is excluded in our model that works in the quenched approximation. The \kx\- 
factorization, known in perturbative physics, has been found to be valid also for the 
non-perturbative dipole-hadron interaction. 

The most outstanding result of this work is the following feature of the QCD 
string that confines the quark and antiquark in the dipole: The QCD string of length 
\f D \ can be exactly represented as an integral over stringless dipoles of sizes £\r D \ with 
< £ < 1 and dipole number density n(£) = l/£ 2 . A similar behavior has been 
observed for the perturbatively computed wave function of a qq onium state in the 
large - N c limit where the numerous emitted gluons inside the onium state represent 
dipoles [36,37]. The decomposition of the string into stringless dipoles has allowed 
us to rewrite the string-hadron scattering process as an incoherent superposition of 
dipole-hadron scattering processes and to extract the unintegrated gluon distribution 
of hadrons and photons, respectively, from our dipole-hadron and dipole-photon 
cross section via \k±\ -factorization. 

We have shown explicitly the microscopic structure of the unintegrated gluon 
distribution in hadrons and photons. For small momenta |fcj_|, the unintegrated 
gluon distributions of hadrons Th{x, k±) are dominated by non-perturbative physics 
and behave as S^/\k±\ where Sh denotes the hadron extension. The l/|fcj_| - behavior 
is a string manifestation as it originates from the linear increase of the total dipole- 
hadron cross section at large dipole sizes. For large momenta, \k±\ > 1 GeV, the un- 
integrated gluon distributions of the hadrons are dominated by perturbative physics 
and show the l/kf_ - behavior induced by the gluon propagator. In the perturbative 
region of large momenta the valence constituents are resolved and the dependence 
of the unintegrated gluon distribution on the hadron extension Sh vanishes. In con- 
trast, the unintegrated gluon distribution of photons depends on the "photon size" 
for large \k^\ because of the Q 2 -dependent normalization of photon wave functions. 

The x- dependence of the unintegrated gluon distribution of hadrons and pho- 
tons has been introduced phenomenologically into our model. Motivated by experi- 
mental observations, we have given a strong energy dependence to the perturbative 
contribution and a weak one to the non-perturbative contribution. Consequently, 
with decreasing x the perturbative contribution increases much stronger than the 
non-perturbative contribution and extends into the small- \k±\ region. A similar 
hard-to-soft diffusion is seen also in the approach of Ivanov and Nikolaev [85] while a 
soft-to-hard diffusion is obtained in the approach of the color glass condensate [86]. 
Considering the integrated gluon distribution of the proton xG p (x, Q 2 ), our non- 
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perturbative contribution dominates for x > 10~ 3 while the perturbative contribu- 
tion becomes relevant for x < 10~ 3 and generates the steep increase of xG p (x, Q 2 ) 
with decreasing x at fixed Q 2 . Also the rise of xG p (x,Q 2 ) with increasing Q 2 at 
fixed x results from the strong energy dependence of the perturbative contribution. 

We have compared the unintegrated gluon distribution of the proton F p {x, k±) 
extracted from our loop-loop correlation model (LLCM) with those obtained from 
the saturation model of Golec-Biernat and Wiisthoff (GBW) [72], the derivative 
of the Gliick, Reya, and Vogt (GRV) parametrization of xG p (x,Q 2 ) [87], and the 
approach of Ivanov and Nikolaev (IN) [85]. For k 2 — > 0, the unintegrated gluon 
distribution of GBW decreases as k 2 and the one of IN as kj_ in contrast to the 
1/ a/AJ - decrease found in our model. In the perturbative region, the LLCM, GRV, 
and IN unintegrated gluon distributions become parallel for x < 10~ 2 and drop as 
l/k]_ with increasing k]_. This perturbative QCD behavior is not reproduced by the 
GBW unintegrated gluon distribution which decreases exponentially with increasing 
k\. The x-dependence of the considered unintegrated gluon distributions is weak in 
the non-perturbative region and becomes stronger as k\ increases. 

We have studied saturation effects of the scattering amplitude in impact param- 
eter space. The impact parameter profiles have been found to respect the black disc 
limit which is imposed by the S'-matrix unitarity. We have computed the impact pa- 
rameter profiles for hadron-hadron and longitudinal photon-proton scattering. They 
have shown the following behavior for the evolution of the opacity and size of the 
interaction particles with increasing energy: The particles become larger and blacker 
as the energy increases. At energies beyond the black disc limit, the opacity satu- 
rates while the expansion of the scattered particles continues. The absolut value of 
the black disc limit and the energy at which it is reached depend on the normaliza- 
tion of the wave functions and the size of the interacting particles, respectively. For 
asymptotic energies a universal result is obtained: the size of each particle increases 
logarithmically with energy as required to guarantee the Froissart bound [7]. 

We have estimated the impact parameter dependent gluon distribution of the 
proton xG(x,Q 2 , \b±\). We have shown a low-x saturation of xG(x,Q 2 , \b±\) as a 
manifestation of the S'-matrix unitarity realized by multiple gluonic interactions. 
The gluon density is found to decrease from the center towards the periphery of 
the proton. The saturation value and the increase of xG(x,Q 2 , \b±\ = 0) towards 
small x depend on the photon virtuality Q 2 . In contrast, at fixed Q 2 , the integrated 
gluon distribution xG(x, Q 2 ) does not saturate because of the growing proton radius 
with decreasing x. In agreement with other investigations of gluon saturation [108, 
111, 112], we have found a slow down of xG(x, Q 2 ) from a powerlike to a squared 
logarithmic rise for Bjorken x values beyond the black disc limit in xG(x, Q 2 , \b±\). 
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The energies at which the saturation effects set in are model dependent. In our 
approach the impact parameter profiles saturate the black disc limit at zero impact 
parameter for y/s > 10 6 GeV in proton-proton scattering and for y/s > 10 7 GeV 
in longitudinal photon-proton scattering when Q 2 > 1 GeV 2 . The saturation of 
xG(x, Q 2 , \b±\) occurs in our approach for Q 2 > 1 GeV 2 at values of x < 10~ 10 
which is far below the HERA and THERA range. Munier et al. [114] have also 
claimed that the black disc limit is not reached in the energy range of HERA. 

We have computed total cross sections a tot for proton-proton, pion-proton, kaon- 
proton, photon-proton, and photon-photon reactions. A stronger energy rise of total 
cross sections has been obtained going from large to small-size particles involved in 
the interaction. This comes, of course, from the strong energy dependence of the 
perturbative component that becomes increasingly important with decreasing par- 
ticle sizes. The growth of the total cross sections becomes weaker for cm. energies 
y/s > 10 6 GeV due to the onset of the black disc limit at = in the profile 
functions. In fact, a transition from a powerlike to an In 2 (s) -increase of total cross 
sections of hadrons and photons is found to set in at y/s ~ 10 6 GeV. For asymptotic 
energies, the total hadronic cross sections become universal, i.e., independent of the 
hadron species considered, and increase in agreement with the Froissart bound [7]. 
We predict for proton-proton scattering at LHC (y/s = 14 TeV) a total cross section 
of <r*°* = 114.2 mb which is in good agreement with the cosmic ray data. 

For differential elastic cross sections da el /dt of proton-proton, pion-proton and 
kaon-proton scattering, the diffraction pattern and also the shrinkage of the diffrac- 
tion peak with increasing energy has been shown in good agreement with experi- 
mental data at small momentum transfers \t\. Around the dip region, where a real 
part is expected to be important, deviations from the data have reflected that our 
T-matrix is purely imaginary. The smaller size of the pion and kaon as compared to 
the proton has become visible in the shift of the dip towards larger values of \t\ in 
the differential elastic cross sections. A differential elastic cross section with a nega- 
tive curvature, C < 0, and a dip at \t\ ~ 0.35 GeV 2 is predicted for LHC. The total 
elastic cross section, obtained by integrating the differential elastic cross section over 
the momentum transfer \t\, has shown similar features as the total hadronic cross 
section discussed above. It takes a value of a el ~ 30 mb at the LHC energy. 

We have also studied the slope parameter B(s) of the forward differential elas- 
tic cross section for proton-proton scattering and the ratios a el /a tot and a tot / B for 
proton-proton, pion-proton and kaon-proton reactions. The slope parameter and 
the ratios are a measure of the squared hadron extension and hadron opacity, re- 
spectively. We have found an In 2 (s) -increase for the slope parameter at large cm. 
energies. The ratios have shown a universal behavior already for cm. energies 
y/s > 10 6 GeV. This behavior reflects most explicitly the saturation of the impact 
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parameter profiles at the black disc limit. The slope parameter at LHC is predicted 
to be B = 21.26 GeV -2 . 

The mentioned observables, namely, total cross sections, the structure function 
of the proton, slope parameters, differential elastic cross sections, elastic cross sec- 
tions and the ratios a el /a tot and a tot /B for proton-proton, pion-proton, kaon-proton, 
photon-proton, and photon-photon reactions are in good agreement with experimen- 
tal data over a wide cm. energy range. We have provided a unified description of 
these reactions using a universal mechanism for the interaction and energy depen- 
dence. This mechanism seems to be appropriate and required by experimental data. 

In spite of the success of our combined perturbative and non-perturbative QCD 
approach to high-energy scattering, there are still several shortcomings and limita- 
tions of our model in its present form that leave room for future improvements: 

• Our approach does not generate energy dependent cross sections because of the 
missed gluon radiation. Thus, for a fundamental understanding of the energy 
dependence of high-energy reactions one has to implement quantum evolution 
in our model analogous to complementary approaches [75]. Attempts aiming 
at a description of high-energy scattering from the QCD Lagrangian have been 
recently reported: In [148] structure functions of deep inelastic scattering at 
small Bjorken x are related to an effective Euclidean field theory. Here one 
hopes that the limit x — > corresponds to critical behavior in the effective the- 
ory. In another recent formalism, the energy dependence of the proton struc- 
ture function has been related successfully to critical properties of an effective 
near light-cone Hamiltonian in a non-perturbative lattice approach [149]. 

• Our model works in the quenched approximation since dynamical fermion 
production is neglected. In addition also gluon vacuum polarization is missed. 
These shortcomings, of course, exclude string breaking in color-dipoles in the 
fundamental and adjoint representation, respectively. To be consistent with 
lattice QCD simulations [26, 80], dynamical fermion and gluon production have 
to be introduced. Suggestions along these lines can be found in [11]. 

• The Gaussian approximation of the functional integrals leads to a dependence 
of cross sections on the choice of the surface that is bounded the Wegner- 
Wislon loops. It is a challenge to find the right surface that is favored by the 
Gaussian approximation. There are several indications supporting the minimal 
surfaces used in this work. Nevertheless, the minimal surfaces are not choosen 
by a variational principle and are not dynamical. Interpreting the surfaces as 
the worldsheets of the QCD strings, this means that our model can neither 
describe quantum fluctuations or excitations of the string nor string flips be- 
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tween two interacting color-dipoles. Accordingly, our model cannot reproduce 
the Liischer term (exitation of the string) in the static quark-antiquark po- 
tential recently confirmed on the lattice [68]. Interesting new developments 
towards a dynamical surface choice and a theory for the dynamics of the con- 
fining strings can be found in [54]. 

• Our ansatz for the gluon field strength correlator does not include explicitly 
the dependence on the path connecting the two gluon field strengths which has 
been recently measured on the lattice [55]. In addition, the non-perturbative 
correlator, being appropriate to describe the scattering between large-size par- 
ticles, overestimates the scattering between small-size particles as shown in 
deep inelastic scattering at high photon virtualities. Moreover, our decompo- 
sition of the correlator into a perturbative plus a non-perturbative component 
is not mandatory. Here it is not clear how to reconcile non-perturbative cor- 
relations with perturbative gluon exchange. Methods complementary to our 
simple two component ansatz can be found in [49,52-54]. 

• The static limit of our model becomes most explicit in the computation of 
the QCD van der Waals interaction between two static color-dipoles. A wrong 
result for the van der Waals potential is obtained from our model because of the 
energy degeneracy between the intermediate octet states and the intial (final) 
singlet states in the static limit. For a meaningful investigation of QCD van 
der Waals forces within our model, one has to go beyond the static limit and 
to describe the limited lifetime of the intermediate octet states appropriately. 

Several other applications and extensions of the model can be investigated in the 
future. The two most natural continuations of this work would be 

• to compute the vector meson production. Some work in this direction has 
already been done. 

• to generalize our model from two to many particle interactions, i.e. from 
the present particle-particle to particle-nucleus or nucleus-nucleus interactions. 
This would mean an extension of the model from the loop-loop correlation 
(WrJC^WrjCy) to multiple loop correlations (WrjCi] ■ ■ ■ W rn [C n }). There 
are already some first attempts in this direction. 
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Appendix A 
Conventions 



For completeness the notational conventions used in this work are listed below. 



A.l Units 

We work in natural units, 

h = c=l. (A.l.l) 
The conversion constant reads he = 197.32705 MeVfm. 



A. 2 Lorentz Vectors 



Following the notations of the textbooks [150,151], we write the contravariant po- 
sition four-vector x^ in its instant form 



x M = x 1 , x 2 , x 3 ) = (t, x, y, z) = (x°, x_l, x 3 ) = (x°, x) 
The covariant position four-vector 

= g^x u = (x , xx,x 2 , x 3 ) = (t, -x, -y, -z) 
is obtained from the contravariant vector by the metric tensor 

f+1 \ 



9tw = 9" 



0-100 
0-10 
\ -1 / 



(A.2.1) 
(A.2.2) 



(A.2.3) 
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Of course, implicit summation over repeated Lorentz indices is understood. The 
scalar product of x^ with four-momentum p M = (p° , p 1 , p 2 , p 3 ) = (E,p) is 



x ■ p = x^p^ = x°po + x pi + x 2 p2 + x 3 p3 = tE — x ■ p . (A. 2. 4) 



A. 3 Light-Cone Coordinates 

For the light-cone coordinates, we employ the convention that has been used by 
Lepage and Brodsky [152] and write the contravariant position four- vector as 



[ X j X j X i J 



(A.3.1) 



with time-like ('light-cone time') and space-like ('light-cone position') components 

-° 1 x 3 and x~ = x° -x 3 , (A.3.2) 



x 



x 



respectively. The covariant position four-vector is again obtained by lowering the 
indices with the metric tensor, x^ = g liu x u , that has the following form in the 
Lepage-Brodsky light-cone coordinates 



/0 2 \ 

2 

0-10 

\0 -1 / 



and 



f° 


1 

2 





\ 


1 

2 

















-1 
















(A.3.3) 



This form was required by the Lorentz-invariance of the scalar product. The scalar 
product of the position four-vector with the momentum four-vector is written as 



x ■ p = x^p^ = x + p + + x p~ + x 1 pi + x 2 p2 = ~(x + p + x p + ) — xj_p± . (A. 3. 4) 

The measure for a four-dimensional space-time integration, d 4 x, expressed in light- 
cone coordinates reads 



d 4 x = dx°dx 1 dx 2 dx 3 = -dx + dx d 2 x±. 

2 



(A.3.5) 
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Appendix B 
Wave Functions 



The light-cone wave functions ipi(zi,fi) provide the distribution of transverse size 
and orientation ri and longitudinal quark momentum fraction Zi to the light-like 
Wegner- Wilson loops W^[Cj] that represent the scattering color-dipoles. In this way, 
they specify the projectiles as mesons, baryons described as quark-diquark systems, 
or photons that fluctuate into a quark-antiquark pair before the interaction. 

B.l Hadron Wave Function 

In this work mesons and baryons are assumed to have a quark-antiquark and quark- 
diquark valence structure, respectively. As quark-diquark systems are equivalent 
to quark-antiquark systems [153], this allows us to model not only mesons but also 
baryons as color-dipoles represented by Wegner- Wilson loops. We use for the hadron 
wave function the phenomenological Gaussian Wirbel-Stech-Bauer ansatz [44] 



1ph{Zi,fi) 



Zj(l ~ Zi) 

2nS 2 h N h 



e -( 2l -i) 2 /(4A 2 2) e -\n\ 2 /(4S 2 h ) 



(B.1.1) 



where the hadron wave function normalization to unity 





requires the normalization constant 




(B.1.3) 
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Table B.l: Hadron Parameters 



Hadron 


Az h 


Sh [fm] 


P,P 


0.3 


0.86 


n ± 


2 


0.607 


K± 


0.57 


0.55 



The different hadrons considered — protons, pions, and kaons — are specified by 
Azh and Sh providing the width for the distributions of the longitudinal momentum 
fraction carried by the quark Zi and transverse spatial extension \fi\, respectively. 
In this work the extension parameter Sh is a fit parameter that should resemble 
approximately the electromagnetic radius of the corresponding hadron [23], while 
Azh = w/(\^2m h ) [44] is fixed by the hadron mass m h and the value w = 0.35 — 
0.5 GeV extracted from experimental data. We find for (anti-)protons Az p = 0.3 
and S p = 0.86 fm, for pions Az n = 2 and S n = 0.607 fm, and for kaons Azk = 0.57 
and Sjc = 0.55 fm which are the values used in the main text. For convenience they 
are summarized in Table B.l. 

Concerning the quark-diquark structure of the baryons, the more conventional 
three-quark structure of a baryon would complicate the model significantly but 
would lead to similar predictions once the model parameters are readjusted [13]. In 
fact, there are also physical arguments that favor the quark-diquark structure of the 
baryon such as the 51 = 1/2 enhancement in semi-leptonic decays of baryons [153] 
and the strong attraction in the scalar diquark channel in the instanton vacuum [154]. 



B.2 Photon Wave Function 

The photon wave function ip 1 (zi,fi,Q 2 ) describes the fluctuation of a photon with 
virtuality Q 2 into a quark-antiquark pair with longitudinal quark momentum frac- 
tion Zi and spatial transverse size and orientation r*. The computation of the corre- 
sponding transition amplitude {qq(zi,ri)\ r y*(Q 2 )) can be performed conveniently in 
light-cone perturbation theory [155] and leads to the following squared wave func- 
tions for transverse (T) and longitudinally (L) polarized photons [69] 

\^( Zi ,n,Q 2 )\ 2 = |^^eH^ + (l-^ 2 ]4^i( e /l^l) + m /^o 2 (^/kll)} (B-2.1) 
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|^(^r1,g 2 )| 2 = |^^eJ{4g 2 ^(l-^) 2 ^o 2 (^/lnl)}, (B.2.2) 

/ 

where a is the fine-structure constant, ej is the electric charge of the quark with 
flavor /, and Kq and K\ are the modified Bessel functions (McDonald functions). 
In the above expressions, 

controlls the transverse size (-distribution) of the emerging dipole, |r*j| oc 1/e/, that 
depends on the quark flavor through the current quark mass nif. 



e 2 f = z i {l-z i )Q 2 +m 2 f (B.2.3) 



For small Q 2 , the perturbatively derived wave functions, (B.2.1) and (B.2.2), 
are not appropriate since the resulting large color- dip oles, i.e. |r*j| oc l/mj ^> 1 fm, 
should encounter non-perturbative effects such as confinement and chiral symmetry 
breaking. To take these effects into account the vector meson dominance (VMD) 
model [156] is usually used. However, the transition from the "partonic" behavior at 
large Q 2 to the "hadronic" one at small Q 2 can be modelled as well by introducing 
g 2 -dependent quark masses, rrif — m,f(Q 2 ), that interpolate between the current 
quarks at large Q 2 and the constituent quarks at small Q 2 [20]. Following this 
approach, we use (B.2.1) and (B.2.2) also in the low-g 2 region but with the running 
quark masses 



m u , d (Q 2 ) = 0.178 GeV ( 1-^-j 6(g 2 >d -g 2 ) , (B.2.4) 

m s (Q 2 ) = 0.121 GeV + 0.129 GeV (^1- ^ 6(g 2 -g 2 ) , (B.2.5) 
and the fixed charm quark mass 

m c = 1.25 GeV , (B.2.6) 

where the parameters Q 2 u d = 1.05 GeV 2 and Q 2 = 1.6 GeV 2 are taken directly 
from [20] while we reduced the values for the constituent quark masses rrif(Q 2 = 0) 
of [20] by about 20%. The smaller constituent quark masses are necessary in order 
to reproduce the total cross sections for 7*p and 7*7* reactions at low Q 2 . Similar 
running quark masses are obtained in a QCD-motivated model of the spontaneous 
chiral symmetry breaking in the instanton vacuum [157] that improve the description 
of 7*p scattering at low Q 2 [158]. 
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Appendix C 
Correlation Functions 



In this appendix we describe explicitly the way from the simple exponential correla- 
tion functions in Euclidean space-time (2.4.10) to their transverse Fourier transforms 
in Minkowski space-time, (2.5.11) and (2.5.16). The first step in this procedure is 
the Fourier transformation of the exponential correlation functions (2.4.10) in four- 
dimensional Euclidean space-time 

D E (K 2 ) = D E (K 2 ) = Jd A ZD E (Z 2 /a 2 )e lKZ 

POO PIT PIT p27T 

= / d\Z\ \Z\ 3 / d(f) 3 sin 2 03 / dfo sin0 2 / dfa D E (Z 2 /a 2 ) e -*l#H3|«**a 
Jo Jo Jo Jo 

a 2 roc -i o_2 

= ^- / d\Z\\Z\ 2 D E {Z 2 /a 2 )J l {\K\\Z\) = — - , (C.0.1) 

\K\ Jo a(K 2 + a- 2 )l 

where J\ is the 1st order Bessel function of the first kind. Here the Euclidean metric 
— Sfj, u and four-dimensional polar coordinates and the corresponding four- volume 
element d 4 Z = d\Z\ \Z\ 3 d(f>3 sin 2 3 (i</>2 sin</> 2 d(f>i have been used. With (C.0.1), one 
obtains 

Now, (C.0.1) and (CO. 2) are analytically continued to Minkowski space-time, K± — > 
ik° or equivalently —5^ v — > g^ v = diag(l, —1, —1, —1), 

D{k 2 ) = , D[(k 2 ) = . (C.0.3) 

a {— k 2 + a~ 2 ) 2 a(— k 2 + a~ 2 ) 2 
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Setting k° = k 3 = 0, which is enforced in the computation of x by 5-functions, one 
finds k 2 = —k\ and consequently 

a(kj_ + a~^)2 a(kj_ + a~ z )2 

The transverse Fourier transformation (2.5.9) of these two expressions is the remain- 
ing step that leads directly to (2.5.11) and (2.5.16). 
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Appendix D 

Non-Forward X-Matrix Element 



In this appendix we calculate the perturbative and non-confining contribution to the 
non-forward (t 7^ 0) T-matrix element. We show explicitly that the non-forward T- 
matrix element depends on the parameters which control the Zi - distribution of the 
wave functions. These parameter are essential for a good description of differential 
elastic cross sections and the slope parameter as shown in Chapter 5. 

The confining contribution to the non-forward (t 7^ 0) T-matrix element is not 
presented. It is much more complicated since some of the integrations cannot be 
performed analytically. Nevertheless, we find numerically that it shows the same 
features concerning the Zj - distributions as the perturbative and non-confining con- 
tributions. 

The perturbative contribution to the non-forward T-matrix element in the small- 
X limit (2.6.3), 

T p (s ,t) = ^ Jd 2 b^ Jdz 1 d 2 r 1 Jdz 2 d 2 r 2 \M^n)\ 2 mz 2 ,r 2 )\ 2 ( X p f , 

(D.0.1) 

reduces upon integration over the impact parameter \b±\ to 
T p (s ,t) = ™p.J d 2 k ± a s (kl) iDfXkl) a,((4 + gl) 2 ) W'® ((k ± + q ± f 

x J dz x j dz 2 H x (zlql) - H x {[z x q\ + 4) 2 ) H 2 (z 2 2 ql) - H 2 {[z 2 q\_ + k A 

(D.0.2) 

with 

Hi {{z t q ± + 4)) := J d 2 r, \^{z h r*)j 2 e** . (D.0.3) 
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and \ipi(zi, fi)\ 2 denoting hadron or photon wave functions. 

The non-confining contribution to the non-forward T-matrix element in the 
small- x limit (2.6.3), 

TZTM = Jd% x e^ jd Zl d 2 ri jdz 2 d 2 r 2 \M^n)\ 2 |^ 2 ,r 2 )| 2 (xD 2 , 

(D.0.4) 

becomes analogously 

rpNp / .x 8zs f tt 2 G 2 (1 - k)\ 2 f d 2 k ± .fs/(2)/L2N .fs/(2) //r . - x 2 \ mn^ 
T „ c (so,*) = -g- ) J (2^p 1 " ^ 1 + J ( D - - 5 ) 

x j\ Zl £dz2 [H x {zlql) - Hi + 4) 2 )] [#2 (^) - #2 ((*2& + 4) 2 ) 

with i7 12 defined in (D.0.3). 

For i = —q\ 7^ 0, both contributions (D.0.2) and (D.0.5) depend on the shape 
of the ^-distribution of the wave function, i.e., for the Gaussian hadron wave func- 
tion (B.l.l), the contributions depend on the width Azh- This Az^-dependence is 
transferred to the differential elastic cross section da el /dt(s,t) given in (5.4.1) and 
its local slope B(s,t) given in (5.3.1). At t — 0, the dependence on the shape of the 
^-distribution of the wave functions disappears because of the normalization of the 
^-distribution as can be seen immediately from (D.0.3). Therefore, in our model the 
total cross sections - related via the optical theorem to the forward (t = 0) T-matrix 
element - do not depend on the parameter that characterize the Zj-distribution of 
the wave function. 
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